SPECTRAL ANALYSIS OF SUBORDINATE BROWNIAN 
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Abstract. We study one-dimensional Levy processes with Levy- 
Khintchine exponent ipi^'^), where i/j is a, complete Bernstein function. 
These processes are subordinate Brownian motions corresponding to 
subordinators, whose Levy measure has completely monotone density; 
or, equivalently, symmetric Levy processes whose Levy measure has com- 
pletely monotone density on (0, oo). Examples include symmetric stable 
processes and relativistic processes. The main result is a formula for the 
generalized eigenfunctions of transition operators of the process killed 
after exiting the half-line. A generalized eigenfunction expansion of the 
transition operators is derived. As an application, a formula for the 
distribution of the first passage time (or the supremum functional) is 
obtained. 



1. Introduction and statement of main results 

In a recent paper |l9], spectral problem for the one-dimensional Cauchy 
process (that is, the symmetric 1-stable process) killed upon exiting a half- 
line or an interval was studied. For the half-line, an explicit formula for 
generalized eigenfunctions of transition operators was obtained using meth- 
ods developed in the theory of linear water waves. The argument of |19] 
relies on some properties specific to the Cauchy process, and it does not 
easily generalize to other Levy processes. The purpose of this article is to 
derive a similar formula for generalized eigenfunctions in a more general 
setting, using a modified method. 

The class of processes considered here consists of symmetric (one-dimen- 
sional) Levy processes Xf with Levy measure having completely monotone 
density function on (0, oo). Alternatively, this class can be described as 
(one-dimensional) subordinate Brownian motions, corresponding to sub- 
ordinators with Levy measure having completely monotone density func- 
tion. Yet another characterization is given by the condition that the Levy- 
Khintchine exponent of Xt has the form ipiC,"^) for a complete Bernstein func- 



tion The equivalence of the above is given in Proposition 2.14 This 



class of Levy processes have attracted much attention in the last decade; 
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see, for example, [Ini[l6l[171|35lll0lllllll3lllll|68l|69]for some recent 
developments. There is also an extensive literature focused specifically on 
symmetric a-stable processes ([2l[3l[9l[IIl[l5l[l9l[20linil36lll7lll8l|5^ 
and relativistic a-stable processes ([121 [131 El EHl |63]), which are included 
here as examples. 

We consider the process Xt killed upon leaving the half-line (0, oo). Let 
p(o,oo) ^ ^(o,oo) and V{A(o,oo)] L'^) denote the corresponding transition op- 
erators, the L°°((0, cxd)) generator, and the domain of .4(o,oo)5 respectively 
(formal definitions are given in the Preliminaries). The following general- 
ization of Theorem 2 from |19] is the main result of the article. Its main 
advantage is the explicit description of the eigenf unctions. 

Theorem 1.1. Suppose that the Levy-Khintchine exponent of Xt has the 
form ipi^"^) for a complete Bernstein function ip. For all A > 0, there is 
a bounded function Fx on (0, oo) which is the eigenfunction of pI^^'°°^ and 

'°°)Fa(x) = e-*^(^')FA(x) and Ao,oo)Fa(x) = -^{X')Fx{x) 
for all X > 0. The function Fx is characterized by its Laplace transform: 

(1.1) CFxiO = , exp - / log—— ^T^rfC 



A2 + e \^Jo e + e "V(A2)-^(e) 
for ^ G C such that Re^ > 0. Furthermore, for x > we have 

(1.2) Fx{x) = sm{Xx + ^x) -Gx{x), 

where -dx G [0,7r/2), andGx{x) is a bounded, completely monotone function 
on (0, oo). More precisely, we have 

n« ,1 r ^ , #(A')(A'-C') 

^' = "1 - 4,(0 

and Gx is the Laplace transform of a finite measure 7a on (0, oo). When 
extends to a function holomorphic in the upper complex half- 

plane {.^ G C : Im^ > 0} and continuous in E C : ImC, > 0}, and 
furthermore 7^ '^/'(A) for all ^ > 0, then the measure 7a is absolutely 

continuous, and 



..4. vrV ^(A2)-^+(-e) 

for ^ > 0. 

In some cases, the above formulae can be substantially simplified. For 
example, when Xt is the symmetric a-stable process, then "i^A = (2 — a)7r/8, 
and Fx{x) = Fi(Ax); see Section |6} 

Any complete Bernstein function extends holomorphically to C \ 

(—00,0]. Hence, the assumption in the last part of Theorem |1.1| concerns 
the existence of boundary values of ip (approached from the upper complex 
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half-plane) along the branch cut (—00, 0]. A formula similar to (1.4) can be 



given in the general discussed in Remark 14.13 



It should be emphasized that although Theorem LI extends the result 
of |l9], its proof is essentially different. 



Remark 1.2. Entirely analytical formulation of Theorem LI is available; 
see [37j for the analytical definition of ^(0,00) • Also, our arguments are 



purely analytical, with the exception of probabilistic proofs of Lemmas 2.10 



and|23Il □ 

The eigenfunctions Fx never belong to L^((0, 00)): this reflects the fact 
that the spectrum of transition operators pI:'^'°°\ considered as operators 
on L^((0,oo)), is purely continuous. Nevertheless, the functions Fx yield a 
generalized eigenfunction expansion of p^'^'^\ In classical eigenfunction ex- 
pansion, a function is decomposed with respect to a complete orthonormal 
set of eigenfunctions, and the operator acts at each component indepen- 
dently. Informally, in generalized eigenfunction expansions, the principle is 
the same, but eigenfunctions are no longer square-integrable and there are 
uncountably many of them. A formal statement is given in the following re- 
sult, which generalizes Theorem 3 in |19], and explicitly provides functional 



calculus for ^(o,oo)- Here .4(o,oo) and V{A[q,oo)] L"^) denote the L^((0, cxd)) 
generator of the semigroup pj^^'°^\ and the domain of ^(0,00) ? respectively 
(again, see Preliminaries for formal deflnitions). 

Theorem 1.3. With the notation of Theorem \l . 1\ let 

POO 

(1.5) n/(A)=/ fix)Fxix)dx for \>0, f ECHO, 00)). 
Jo 

Then ^j2fnJ\ extends to a unitary operator on L^((0, 00)), and 
npf ■°^V(A) = e-*^(^')n/(A) for f G ^^((0, 00)). 

Furthermore, f G 1^(^(0,00) ; -^^) if and only if ilj{X'^)Ilf{X) is in L^((0, 00)), 
and 



nAo,oo)/(A) = -ij{X')Uf{X) for f G I?(Ao,oo); L 



2^ 



Remark 1.4. Below we provide a complete proof of Theorem 1.3 under 
an additional condition that 11 is injective, and we verify that 11 is in- 
jective when Xt is the symmetric a-stable process, a G (0,2] (that is, 
'^{0 = ^°''^)- In the preliminary version of this article, a relatively easy 
to check, but rather restrictive sufficient condition for injectivity of 11 was 
given (Lemma 19 in [S3]), and it was conjectured that in fact 11 is always 
injective. After a few months, this conjecture was solved in affirmative 
by Jacek Malecki, and the proof of injectivity of 11 is given in the recent 
preprint |56]. We announce the result of [56] and state Theorem 1.3 in full 
generality. □ 



Remark 1.5. There is a striking similarity between formulae (1.1 ) and (1.4) 
on one hand, and some formulae in fluctuation theory of Levy processes on 
the other. More precisely, the exponent in (1.1) and (1.4) resembles known 
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formulae for the Laplace exponent of ladder processes (see Corollary 9.7 
in [28j) and the Laplace transform of the supremum functional (Theorem 1 
in [4j). This phenomenon can be easily explained at the analytical level: in 
both cases explicit formulae are obtained using the Wiener-Hopf method, 
a technique for solving some integral equations in the half-line via Fourier 
transform. Fluctuation theory of Levy processes relies upon the Wiener- 
Hopf factorization of the A-potential operator (that is, the resolvent) of the 
transition semigroup of a Levy process. The Fourier symbol of this operator 
is 1/(A + ipl^"^)). In the proof of Theorem the operator with Fourier 
symbol {ip{X'^) — ^('C^))/(A^ — is factorized. 

A probabilistic formulation of the fluctuation theory is available through 
local times, excursion theory and ladder processes; see jH [231 [571 |6l] ■ It is an 
interesting open problem whether there is a similar probabilistic derivation 
of the formula for Fx. □ 



The proof of Theorem L3 (under the assumption that LI is injective) 
follows the line of |19]. Proofs of Theorems |L 1 1 and 1.3 are rather technical. 



Therefore, it may be helpful to keep in mind the following well-recognized 
example. 

Example 1.6. Suppose that Xt is the Brownian motion, with variance 2t 
(so that the generator of Xt is the one-dimensional Laplace operator (P /dx^, 
and ip{C) = 0- I^ this case ^(o,oo) is the Laplace operator in (0,oo) with 
Dirichlet boundary condition at 0, and p^^'^'°°-' is the classical heat semigroup 
on (0, cxd) with Dirichlet boundary condition at 0. The eigenfunctions of 
^(o,oo) and pI^^'°°^ are simply Fx{x) = sin(Ax), with corresponding eigenval- 
ues — A^ and e~*^^, respectively. The integral transform n is the Fourier sine 



transform on (0, oo), and Theorem 1.3 states that the Fourier sine transform 



is the unitary mapping (up to a constant factor A/2/7r), which diagonalizes 
the action of the heat semigroup. 



To the author's knowledge, no results similar to Theorems |1.1| and |1.3| 
were available for Levy processes other than the Brownian motion (possibly 
with drift; this requires a minor modification to the above example) and 
the Cauchy process (studied in [IH])- Although generalized eigenfunction 



expansions similar to Theorem 1.3 have been studied for various classes of 
operators (see [30] for the case of general Markov processes, and [6T1 [67] for 
diffusion processes with Feynman-Kac potential), the explicit description 
of eigenfunctions was not available. Consequently, their applications were 



limited. The novelty of Theorem 1.3 is in that the generalized eigenfunc 



tions Fx{x) are given in a fairly explicit form, allowing various estimates, 
asymptotic analysis and numerical approximation. For this reason, one can 



expect that Theorems 1.1 and 1.3 will find a variety of applications; some 
are already discu sse d be low. 

yield formulae for the transition density pf''°°\x, y) 



Theorems 



1.1 



and 



1.3 



of Xt killed upon leaving the half-line (0, oo), and for the distribution of the 
first exit time T(o,oo) from the positive half-line. These results are contained 
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in the following two theorems. By we denote the probability for the 
process Xt starting at a fixed point x > 0. 



;i.6) 



Pt 



(0,00) 



Theorem 1.7. Let t>0. If e 

2 

7T 

Theorem 1.8. Lett>0. If 



tV'(C^) j^g integrable in ^ > 0, then 

POO 

/ e-'^^^"^Fxix)Fxiy)dX forx,y>0. 
Jo 



(1.7) 

and 

(1.8) 

then 

(1.9) P.(r(o,oo) > t) 



sup 



< 2, 




< 00, 



for X > 0. 



If \/V'(C^)'0'(^^)e zs integrable in > 1 for all t > 0, then furthermore 

(1.10) P.(r(o,oo) e rft) ^ 
/or x,t > 0. 



IT 



^/^'{X^)iJ{X^) e-'^^^'^Fx{x)dX ] dt 



Remark 1.9. Theorem 1.8 in full generality is proved in the recent 
preprint jHS]- In this article we provide a much simpler proof, under a 
more restrictive condition 



(1.11) 



hmsup 



< 1, 



lim sup 



< 1. 



This condition implies both (1.7) and ( 1.8[ ), as well as integr ability of 
■\/ {i'^)e~^^'^^ ) (^gge the proof of the theorem). The proof given be- 
low relies on Theorem 1.3 Noteworthy, many interesting examples, in- 
cluding symmetric a-stable processes and relativistic a-stable processes, 
satisfy (1.11). By Corollary 4.26, condition ( |1.11[ ) is satisfied when if) is 
regularly varying of positive order at and at 00. □ 



It should be emphasized that the integrand in (1.6) has two oscillatory 
factors, and cancellations in the integral are essential. Hence, formula (1.6) 
is problematic for numerical computations. For the Cauchy process, a sig- 
nificant simplification of (1.6) is possible (see Theorem 4 in |:49j). It is an 
open problem whether formula (1.6) can be simplified in the more general 
case, even for symmetric a-stable processes for general a G (0,2]. 

Cancellation is of less importance in Theorem 1.8, see ^56j, where (1.10) 
is used to obtain estimates for the density function of the distribution of 

T(0,oo)- 

Since ^|V'"(OI/V''(0 < 2 for all > 0, the supremum in (1.7) is always 
not greater than 2 (see Preliminaries). The assumptions (1.7) and (1.8) are 
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rather mild regularity and growth conditions for examples are given 

in Section [6l 

Remark 1.10. By symmetry, the first exit time T(o,oo) for Xf starting at 
a; > has the same distribution as the first passage time r^. through a 
barrier at x for Xf starting at 0. On the other hand, first passage times are 
related to the supremum functional: 



'0 



sup < a; = PoiTx > t) = Px(r(o,oo) > t). 
se[o,i] / 



Hence, Theorem L8 gives, in a rather general setting, an explicit expression 
for the distribution of first passage times and the supremum functional, 
which are fundamental objects in fluctuation theory of Levy processes. 
Noteworthy, the double Laplace transform (in t and x) of the distribution 
of Tx is known for general Levy processes since 1957 (Theorem 1 in [4j). 
Nevertheless, explicit formulae for Px(''"(o,oo) > t) are known only in some 
special cases; see |2H l32[ [351 ISTl [52] for some recent developments in this 
area. Also, a formula for the single Laplace transform (in t) of P(r3; > t) for 
a large class of symmetric Levy processes was obtained only very recently 
in [Sg. □ 



Theorems 1.1 and 1.3 can be applied to certain systems of PDEs, related 
to traces of two-dimensional diffusions, as discussed in Section [7j Further- 
more, there are at least four recent preprints exploiting the results of the 
present article. An application of Theorem |L1| to the spectral theory of 
symmetric a-stable process in the interval can be found in [54j . The case of 
the relativistic 1-stable process in half-line and interval is studied in detail 
in |38]. This application may be of particular interest in quantum physics; 



see |26[ |29| [58] for related research. In [27], Theorems |1.1| and 1.3 



are 



used to obtain refined semi-classical asymptotics for eigenvalues of higher- 
dimensional isotropic a-stable processes in domains. Finally, the proof of 



Theorems 1.3 and 1.8 in full generality, as well as their application to a 
detailed analysis of first passage times, can be found in [56j. 

We conclude the introduction with a brief description of the structure 
of the article. The Preliminaries section contains definitions and standard 
properties used in the article, and auxiliary lemmas. It includes some back- 
ground on distribution theory. Levy processes and their generators, subor- 
dinate Brownian motions and complete Bernstein functions. In Section |3] a 
certain transformation related to the Wiener-Hopf factorization is studied. 
Next, in Section [4], we derive the formula for Fx and prove Theorem [1.1[ 



Theorems 1.3, 1.7 and 1.8 are proved in Section |5j Examples, including 



symmetric a-stable processes and relativistic a-stable processes, are stud- 
ied in detail in Section [6j In Section [7] we give an application of our results 



to systems of PDEs (Theorem 7.1) 



2. Preliminaries 



2.1. Distribution theory. In the theory of partial differential equations, 
it is a common case that one first finds a weak solution to a problem, and 
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then, by showing that the solution is sufficiently regular, one argues that it is 



in fact a strong solution. This approach will be used to prove Theorem 1A_ 
first we find a distributional eigenfunction of ^(o,oo), sjid then we show that 
in fact it belongs to the domain of ^(0,00) • Our argument relies heavily 
on Fourier methods, and the operators involved are typically non-local, so 
the required background on distribution theory is slightly different from the 
one used in partial differential equations. It should be emphasized that the 
use of distributions could be avoided, at least partially, at the price of less 
clear exposition. However, the language of distribution theory seems to be 
well-suited for problems involving generators of killed Levy processes. 

Let S denote the class of Schwartz functions in R*^ {d = 1,2,...), and let 
S' be the space od tempered distributions in R'^. If G 5 and F G S', 
we write {F, if) for the value of F at ip. Below we recall some well-known 
properties of tempered distributions; for a detailed exposition of the theory, 
see e.g. [71] . 

If (f E S, the Fourier transform of ip is J^^piC.) = j-^d e^^'^ip{x)dx. For 
F G S', J-'F is the tempered distribution satisfying {J^F,ip) = {F,J^(p). 

The convolution of (pi,(p2 ^ S is defined in the usual way, (pi * <p>2{x) = 
Jiid ^i{y)V2{x — y)dy. When ip) E S and F G S\ then F * is an infinitely 
smooth function, defined by 

F * ip{x) = {F, ip^), where ip^{y) = ip{x - y). 

The convolution of two distributions is not well-defined in general. Suppose 
that Fi,F2 G S'. We say that Fi and F2 are S'-convolvable if for all ipi, ip2 G 
5, the functions Fx * ip\ and F2 * v^2 are convolvable in the usual sense, i.e. 
the integral J-^d{Fi * ipi){y){F2 * ip2){x — y)dy exists for all x. When this is 
the case, the S' -convolution Fi * F2 is the unique distribution F satisfying 
F * [ipi * ip2) = {Fi * ipi) * [F2 * ip2) for ipi,ip2 E S. Note that there are 
other non-equivalent definitions of the convolution of distributions, and S'- 
convolution is often denoted as Fi ® F2; for the discussion of various notions 
of convolvability, the reader is referred to |22| 171] . 

Recall that the support of a distribution F is the smallest closed set 
suppF with the property that {F, ip) = for aA\ ip E S such that ip{x) = 
for X G supp-F. If any of the tempered distributions Fi, F2 has compact 
support, or (in the one-dimensional case) if both Fi and F2 are supported 
in [0, 00), then Fi and F2 are automatically jS'-convolvable. 

It is well known that the distributions {Fi * F2) * F3 and Fi * {F2 * F^) need 
not be equal; however, if the pairs (^1,^2) and (^2,^3) are iS'-convolvable, 
and furthermore the functions Fi* ipi, F2* ip2, and F3 * ip^ are convolvable 
(that is, (Fi * ipi){y){F2 * ip2){z){F3 * ip3){x — y — z) is integrable in y, z for 
all X G R'^) for any y92, 9^3 £ 5, then the ^'-convolution of Fi, F2 and F3 
is associative; see [TT], Section 4.2.8. 

Any iS'-convolvable distributions Fi, F2 satisfy the exchange formula 
J-{Fi * F2) = J^Fi ■ J-'F2, where the multiplication of distributions TFi and 
J^F2 extends standard multiplication of functions in an appropriate manner, 
see [Ml 171] . Since we only use the exchange formula when J^Fi and J-'F2 
are genuine functions, or when J^Fi is a measure and J-'F2 is a function. 
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we do not discuss the notion of multiplication for general distributions and 
refer the interested reader to [39| l65t \7T\ . 

In the one-dimensional case d = 1, the Laplace transform can be defined 
for tempered distributions. If G 5 is supported in [0, oo), then the Laplace 
transform of ip is denoted by = J^e-^''ip{x)dx (Re^ > 0). This 

is a holomorphic function of ^ in the right complex half-plane Re^ > 0, 
continuous at the boundary. Clearly, J^ip{^) = C(p{—i^). If F E S' is 
supported in [0, oo), the Laplace transform of F is defined for ,^ G C with 
Re^ > by 

CF{0 = {F,e^h), 

where e^{x) = e~^^ and h is any infinitely smooth function such that h{x) = 
1 for X > and h{x) = for x < —1. This definition does not depend on 
the choice of h (see formula (9.1.4) in |71)). The above definition of CF 
extends the usual definition of the Laplace transform of signed measures on 
[0,oo). 

Note that while the Fourier transform of a distribution is again a distri- 
bution, the Laplace transform is always a (holomorphic) function. For a 
fixed if: > 0, the function CF{t — is) is the Fourier transform of the distri- 
bution CfF (where again ej(x) = e~*^), and as t — )• O"*", the distributions 
£F{t — is) converge to J-'F in S'. The exchange formula holds also for the 
Laplace transform: if both Fi,F2 G S' are supported in [0, oo), then they 
are 5'-convolvable and £(Fi*F2)(0 = CFi{C)CF2{0 for all ^ with Re^ > 0. 

We say that two distributions Fi, F2 restricted to an open set D are equal 
if {Fi, ip) = {F2, ip) for any ip E S vanishing in Il'^\D. Equivalently: F1 — F2 
is supported in R"' \ D. 

2.2. Transition semigroups and generators. Below we recall and ex- 
tend some standard definitions and properties, which can be found, for 
example, in IlH |5l |6l]- In this article we are only concerned with one- 



dimensional subordinate Brownian motions. However, Lemma 2.10 below 
might be of interest for general Levy processes. For this reason, in this and 
the next subsection, we work with the general case. 

Let Xt be a Levy process in R'^. We write P^^ and E^: for the probability 
and expectation for the process Xt which starts at x G R'^. The process Xt 
is completely determined by its Levy-Khintchine exponent \E': for ^ G R*^ 
we have 

(2.1) Eoe'«-^* = e-**(^\ 
where 

(2.2) M/(0 = /3e ■ e - ^7 ■ e + / (1 - e'^-' + l\,\<ii^ ■ z)u{dz). 

Here /3 is a nonnegative definite d x d matrix (the diffusion coefficient) , 7 
is a vector in R*^ (the drift), and u is the Levy measure of Xf. a Radon 
measure on R'^\ {0} such that j-^^ min(l, \z\'^)i>{dz) < 00. If z/ is absolutely 
continuous, we denote its density function with the same symbol z/(z). 
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The transition operators of Xt are defined by 



PJ{x) = EJiXt 



f{y)P,{Xt e dy), 



t > 0, X eR, 



wfienever tlie integral is absolutely convergent. Each Pf is a convolution 
operator, with convolution kernel given by the distribution of {—Xf) under 
Pq. By the Levy-Khintchine formula (2.1), 



(2.3) Pte^ = e-**(«)ee 
Furthermore, for ip & S, 



for 



(2.4) 



J^PMO = e-**(«)-Fy.(0 



see Theorem 3.3.3 in [T] (note that in [1], Fourier transform is defined with 
Q-^S.-x instead of e*^'^). The operators Pt form a contraction semigroup on 
each of the Lebesgue spaces L'^{Il'^) {p G [l,cx3]), on the space Cb(R^) 
of bounded continuous functions, and on the space Co(R'^) of continuous 
functions vanishing at infinity (all equipped with usual norms). On L^CR'^) 
{p G [l,oo)) and on Co(R''), this semigroup is strongly continuous (see 
Section 3.4 in [Tj). The generator of each of these semigroups is denoted by 
the same symbol A, and we write 'D[A] X) for the corresponding domain, 
where X is either Lp(R^) {p G [l,oo]), Cb(R'^) or Co(R'^). More precisely, 
/ G V{A; f eX and the hmit 



Af 



lim 



Ptf-f 
t 



exists in the topology of X. Note that the limit, if it exists, does not 
depend on the choice of X (apart from the fact that in Cq(R,'^) and Ciy(Rf') 
it is defined pointwise, while in ^^(R'^) only up to a set of zero measure). 
Therefore, using a single symbol A for operators acting on different domains 
VIA] X) causes no confusion. The Schwartz class 5 is a core of A on each 



of the spaces ^^(R'^) {p G [1, oo)) and Co(R'^), and by (2.4), for G 5 



(2.5) 



^^^(0 = -^(O^^(e), 



^ G R. 



We abbreviate, for example, V{A] Lp{R'^)) to V{A] L^). 



Remark 2.1. 



The spectral theory of Pt and A is very simple, thanks to 
the Levy Khintchine formula (2.3). Indeed, the function e*^'^ is the eigen- 
function of Pt, with eigenvalue e"***^^-*. Hence, e*^'^ belongs to T){A; L°°), 
and it is an eigenfunction of A with eigenvalue — \l/(,^). The generalized 
eigenfunction expansion of Pt and A is given by ( 2.4^ and (2.5); the Fourier 
transform plays the same role as the integral transform 11 in Theorem |1.3 

If Xt is symmetric, then ^^(0 = ^(—0 real. Since sin(^ ■ x + i)) is a 
linear combination of e*^'^ and e**-"^^'^, it is also the eigenfunction of Pt and 
A. Note that the eigenfunctions F\ in Theorem lA behave asymptotically 
as sin(Aa; + -^a) as x — )■ oo. □ 
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When ip G S, then the generator of Xt can be written in the form 
A(p{x) = l3V(p{x) ■ V(p{x) + 7 ■ V(p{x) 

(2.6) 



+ / {(p{x + z) - (p{x) - l\z\<iz ■V(p{x))i^{dz), 



where /3, 7 and u are as in (2.2) (Theorem 31.5 in [6 
(2.7) 



L°°(R<*) 



). Furthermore, 



where C depends only on the process Xt, and ||v'||c2(Rd) is the maximum 
of L°°(R'^) norms of if and its first and second partial derivatives (see the 
proof of Theorem 31.5 in 



We need the following extension of (2.6). For Co(R ) instead of C;,(R ), 
this is well-known, see Theorem 31.5 in [64j. The statement given below 
is, however, difficult to find in the literature, so we provide a proof. Let 
C^(R'^) be the class of functions F such that F and all partial derivatives 
of F (of all orders) belong to Cfe(R'^). 

Lemma 2.2. If F e C^(R'^), then F G V{A; Cb) and 
AF{x) = /3VF{x) ■ VF{x) + 7 ■ VF{x) 

(2.8) 



+ 



{F{x + z) -F{x) 



R'' 



-kl<l 



z-VF{x))iy{dz). 



Remark 2.3. Note that in contrast to the Co(R'^) case, F may fail to 
belong to 'DlA; Cb) when F and its first and second partial derivatives are 
in Cb(R.'^). A simple counterexample can be easily constructed for one- 
dimensional Brownian motion: if F, F', F" G Cb(R), but F" is not uniformly 
continuous, then the convergence of {Ptf — f)/t may fail to be uniform. We 
omit the details. □ 



Proof of Lemma 2^. Let ||/||c*^(R'«) denote the maximum of L°°(R ) norms 
of / and its partial derivatives of order not greater than /c, and denote M = 
ll-^llc*(R'*)- There is a sequence (fn^<S such that (fn and its first and second 
partial derivatives converge locally uniformly to F and the corresponding 
first and second partial derivatives of F, and furthermore ||v5n||c*{R'') 

< CM 

(with C depending only on the dimension d). For / G V{A; L°°) and 
x G R*^, we have (see [25j ) 



(2.9) 



Ptfix)-fix) 



PsAf{x)ds 



APsf{x)ds. 



In particular, \\Ptf - fWi^^Rd) < t||^/||ioo(R,d). Since (fn ^ S C V{A;Co) 
by (2.9) we have 



Pt'-Pnix) - ipn{x) 



Aifr 



< 



APs(Pn{x)ds - AiPn{x] 
\A{PsiPn - (Pn)ix)\ds. 
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By (2.7) 



t 



< 



t 



\PsVn - V^n\\c^(R<l)ds. 



Since Pg and A commute with partial differential operators on S, by (2.9) 



t 



Aipnix) 



t 



-J I ^\\AVn\\cl{-R.'i)ds = - \\A(pn\\c^{B,dy 







Finally, by (2.7) 



Pt(Pn(,x) - i^nix) 



- ALPn{x] 



t CtM 

< 2 ll'^"IICb(R'') - 



As n — i- oo, we have ipn{x) — ?■ F{x) and Pt(fn{x) — ?■ PtF{x) (by domi- 
nated convergence). By (2.6), Taylor's theorem and dominated convergence, 
Aipn{x) converges to the right-hand side of (2.8), which we denote by G{x), 
and the convergence is locally uniform in x G R''. It follows that 



PtF{x) - F{x) 
t 



G(x) 



< 



CtM 



Since x G R'^ was arbitrary, we have F G 'D{A; Ch) and AF{x) = G{x), as 



desired. 



□ 



Informally, the generator A is also a convolution operator, but the con- 
volution kernel is a tempered distribution. We give this a precise meaning; 
see also [SlEl], and |8] for the case of symmetric a-stable processes. 

Definition 2.4. The distributional generator of Xt is the tempered distri- 
bution A E S' defined by the formula 



{A,^)=A^{0), 



Let (p{x) = ^p{—x) for E S. We define A (the reflection of A) by the 
formula {A, if) = {A, ip). Hence, {A, if) = A*ip{0) and A(p = A*ip for if E S 
(see m]). By (g, FA{^) = and J^A(0 = -^(0 for ^ G R"^. 

Proposition 2.5. When F G V{A; L°°), then AF = A*F. 

Proof. We have, by the definition of iS'-convolution, 

(A * F) * Lfi * (^2 = (A * (^i) * (F * Lp2) = (A^Pl) * F * (f2- 

Recall that S C V{A]L^). Hence, A(fi is the ^^(R'^) limit of {Ptipi - 
(pi)/t. Since ipi, ip2 and the convolution kernel of Pt are integrable, and F 
is bounded, by Fubini we have 

{Ptipi - ip^) * F * (p2 = {PtF - F)*ipi* ip2- 

By dominated convergence, 

{PtF - F) * v?! * (^2 



(A * F) * 091 * 099 = lim 



{AF) * V9i * 



as desired. 



□ 
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Corollary 2.6. If F e C^(R'^) (as in Lemma \2^, then AF e C^(R'^), 
and A commutes with partial differential operators on C^(R'^). 

Proof. Led D b e a partial derivative operator (of arbitrary order). By 
Proposition we have DAF = D{A* F) = A* (DF) = ADF. □ 

When Xt is one-dimensional and symmetric, then (2.2) takes the form 

/oo 
il-cosi^z))iyidz), eeR, 
-oo 

where /3 > and mm{l, z'^)iy{dz) < oo. For / G C^(R'^) (as in 
Lemma 2.2), we have 

(2.11) Af{x) = f3f"{x) + pv / {fix + z)- f{x))u{dz), xeR. 

Jr. 

Here pv J denotes the Cauchy principal value integral: 



pv / (fix + z)- fix)Mdz) = hm / (fix + z)- f{x)Mdz). 

Jr ^^0+ JR\{-e,e) 

Of course, when Xf is symmetric, then A = A. 

2.3. Killed process and its generator. The main references for the no- 
tion of a killed process (or part of a Markov process) are [7l[25|, where gen- 
eral strong Markov processes are studied. Here we consider a Levy process 
Xf in R'^. Let T£, be the time of first exit from D, Td = inf {t > : Xf ^ D}. 
The killed process X^ is a strong Markov process in D with lifetime r^j, such 
that X^ = Xt for all t < td. The transition operators P/^, the generator 
Ad, and domains T>[Ad]'^) corresponding to the killed process X[' are 
defined in the same way as for the free process Xt, that is. 



P,^/(X) = E./(Xf ) = E,{fiXt)lt<rn), 



and 



Aof= limj^^^^, 

whenever the limit exists in the topology of the function space X; in 
this case, we write / G T>{Ad'-,'^)- Again we abbreviate, for example, 
V{Ad;Lp{D)) to V{Ad; LP). 

When Xt is the Brownian motion, then A is the Laplace operator in R'^, 
and Ad is the Laplace operator in D with Dirichlet boundary condition. 
Even in this case the relation between Ad and A is very delicate. On the 



other hand, this relation is crucial for the proof of Theorem 1.1 Hence, we 
now discuss this topic in more detail. 

Suppose first that Xt is a compound Poisson process: /3 = 0, 7 = and u 
is a finite measure. In this case, the situation is clear. The following result 
seems to be well-known to specialists, but it is difficult to find in literature 
(cf. Example 3.3.7 in [Ij for the case D = R^). 
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Lemma 2.7. Let Xt be a compound Poisson process in R'^, with Levy mea- 
sure V, and let p G [1,cxd]. Then the semigroup Pf is strongly continuous 
on L^{D), the generator Ad is a bounded operator on each L^{D), and 

Ad fix) = lD{x)Afix) = 1d{x) [ {f{x + z) - f{x))v{dz) 

for all f G LP{D), with the usual extension f{x) = for x ^ D. 

Proof. Denote by M = u{Il'^ \ {0}) the intensity of jumps of Xt, and let 
< Ti < T2 < ... be the sequence of (random) times when Xt jumps. Let / 
be in some L^[D) {p G [1, oo]), and let x G -D be the starting point of Xt. 
Note that 

P,^/(x) = E,(/(X,)1,<,J = E.(/(X,)1,<,J 

+ E,{f{Xt)l 

Tl<t<T2 ) 

+ Ea;(/(Xt)1^2<(<^^); 

indeed, in the first summand t < Ti implies that t < Td, and in the second 
one, when ti < t < r2 and t > td, then f{Xt) = 0. We have f{Xt) = f{x) 
in the first summand, and in the second one, f{Xt) = /(-^rj is independent 
of (ri,r2). Furthermore, Px{Xtj^ G x + dz) — M~^v{dz). Hence, if v{E) — 
v{—E), we have 



E, 



JiXr,)^ [ fix + z)u{dz)^ [ f{x-z)i>{dz)^f*i>{x). 



It follows that, 

Pt^'fix) = f(x)F,(t < Ti) 

+ M- V * i>{x)P,{n <t<T2) 

+ ^x{f{Xt)lT2<t<TD)j 

and so, finally, 

Pt^'fix) - fix) 



{f*u{x)-Mf{x)) 



(2.12) < 



Px{t <Ti)-l + Mt 



+ 



Px(Ti <t<T2)-Mt 



+ 



Mt 

^x{f{Xt)lT2<t<TD) 



f * i>{x) 



t 



The proof will be complete once we show that each of the summands on the 
right-hand side converges to in L^{D). 

Since ri,T2, ... are increase times of a Poisson process, we have 



Px(Tfc < t < Tfe+i) 



k\ 



t>0,k^ 0,1, 



where tq = 0. In particular, |P^(t < ri) - 1 + Mt| = je"^^* - 1 + Mt\ < 
{Mt)^/2 and |P^(ri < t < T2)-Mt\ = \Mte-^^-Mt\ < {Mt)'^, which shows 
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that the first two summands in (2.12) converge to in L^CR,'^). Furthermore, 

T2<t<TD 

)|<E,(|/(Xi)|l 

T2<t) 

OO 

= y^^E^df {XrJ\lrk<t<Tk+i) 
k=2 

oo 

= 5^E,|/(X.J|P,.(r, <t<r,+i) 

k=2 

oo 



I/I * {u*^){x) {Mtf 



Mt 



k=2 



k\ 



Since the total mass of i)* is M , we conclude that 



|Ex(|/|(X,)l 



r2<t<TD)\\ LP{D) 



< 



k=2 

oo 



LP(R'*) 



< 



E 

k=2 



k\ 

mf -Mt 

(1 - e-*^* - Mte-''') 



LP^R.'i) 



which is of order as t — )■ 0^. Hence, also the third summand in (2.12) 
converges to 0. □ 

When Xt is not a compound Poisson process, then the generators of Xt 
and X^ are unbounded operators, and the relation between the domains 
of A and is much less obvious. A point x G dD is said to be a regular 
boundary point of D if inf{t > : ^ D} = a.s. with respect to P^. (note 
that here the inequality t > is strict, while t > is used in the definition 
of T£)). Let Cq{D) denote the space of Cq(R'^) functions vanishing in R^\Z}. 
If every x G dD is a regular boundary point and Xt has the strong Feller 
property (that is, Pf maps L°°(R'^) to C;,(R°')), then the operators form 
a strongly continuous contraction semigroup on Co{D) (that is, Xf has the 
Feller property; see [IE])- Note that if Xt is a symmetric Levy process in 
R (and not a compound Poisson process) and D = (0, oo), then 0, the only 
boundary point of D, is regular (see Theorem 47.5 in [64j). 

We have the following fundamental result due to Dynkin. 

Definition 2.8 (Dynkin characteristic operator, fI^). Let x G R'^ and 

denote by r^.^^ the first exit time from the ball centered at x and with radius 
e. If the limit 

EJ{X,)-f{x) 



21./ 



lim 

£-5-0+ 



E r 

' x,e 



exists, we write / G I?(2t.). If / G T'(2t. 
and we define Qlofix) = St^;/. 



for all X G -D, we write / G ©(St/)) 



Lemma 2.9 (Theorem 5.5 in [25] for Levy processes). Suppose that Xt is 
a Levy process in H'^ which is a strong Feller process. Suppose furthermore 
that D C R'^ is open, and every boundary point of D is regular. Let f G 
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Co{D). Then f e V{Ad; Cq) if and only if f e ©(SId) and ^of e Co{D). 
In this case, Aofix) = 2tD/(x) for x G -D. □ 

The above lemma states that, in a sense, the operators Ad and A have 
the same pointwise definition, but their domains are different. This can 
be made even more exphcit using the language of distribution theory, as 
we now d escri be. Recall that A is the distributional generator, defined in 
Definition 2A^ and A is the refiection of A. 

Lemma 2.10. Suppose that Xf is a Levy process in R'^ which is a strong 
Feller process. Suppose furthermore that D C R*^ is open, and every bound- 
ary point of D is regular. Let F G Cq{D). If the distribution A * F 
restricted to D is equal to a Cq{D) function, then F G T>{Ad'-,Cq) and 
AdF{x) = a* F{x) for xe D. 

In the proof, we use the Dynkin's formula ([25j, formula (5.8)): when r 
is a Markov time, E^t < oo, and / G I^IA; Co), then 



(2.13) 



EJ{X,) - fix) 



AfiXs)ds 



X G R"". 



We denote by B{x,6) and B{x,e) the open and closed balls centered at x 
and with radius e. 

Proof of Lemma 2.1C\ Fix x G -D and let r > 0. Let Qr G iS be an approxi- 
mation to identity: Qriy) > 0, j-^d9r{.y)dy = 1 and gr{y) = when \y\ > r. 
We define f = F * g^. Then / G ^^^(R'^) C V{A\ Cb) (see Lemma and 

(2.14) Af{y) =A*{F* (?,) = {A*F)* g,. = F*{A* g^). 

(The convolution of A, F and gr is associative because F is a bounded 
function, g^- has compact support, and A * 99 is integrable for all G S.) In 
particular, Af G Co(R'^), and therefore / G Vi^A] Cq). Let Tx,e be the time 
of first exit from B{x,e). By Dynkin's formula (2.13), 

E./(X(r,,e)) - fix) 



E r 

' x,e 



Afix) 



E r 

' x,e 



E, 



{Af{Xt)-Af{x))dt 







< sup \Af{y)-Af{x)\. 

y&B{x,e) 



Suppose that r + e < dist(x, R'^ \ D). Then using (2.14) for the right-hand 
side, we obtain 



E./(X(r,,,)) - fix) 



Afix) 



< sup 

y&B{x,r+e) 



\A*Fiy) - A*Fix)\. 



Consider the limit as r — > 0^. Then g,. approximates the Dirac delta mea- 
sure Sq. Hence, f = F * g^. converges uniformly to F, and (again by (2.14 
Afix) converges to A* Fix) (x G -D is fixed). It follows that 

E,F(X(r,,,)) - Fix 



Ex'Tx e 



A* Fix) 



< sup 

y&B{x,e) 



\A*Fiy) - A*Fix) 
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As £ — )■ 0^, the right-hand side tends to zero. Therefore, F 6 Vi^x) and 
= A* F{x). Since x G -D was arbitrary, we conclude that F G VjO li)), 
and 2tD-F(a;) = A*F{x) for x E D. The result follows from Lemma 



2.9 



□ 



Lemma [2. 10| seems to be now. However, a related distributional approach 
to generators of the free and killed processes was used for example in [8j, 
in the study of harmonic functions with respect to symmetric a-stable pro- 
cesses. 

We need the following auxiliary result. Its analogue is relatively easy to 
prove using Dirichlet forms. The uniform version seems to be new, especially 
in the case of unbounded D and F not vanishing at infinity. 

Lemma 2.11. Suppose that r > 0, D C R'^ is an open set, F G C^(R'^), 
and F{x) = whenever dist(x,R'^ \ D) < r. Suppose furthermore, that 
AF{x) = for all x e dD. Then F G V{Ad; Cb) and AdF{x) = AF{x) 
for X E D. 

Recall that B{x,e) is the open ball centered at x and with radius e. In 
the proof, the following well-known result is used. For completeness, we 
give a simple proof. 

Proposition 2.12. For all e > 0, Po(t > TB{o,e))/'t is bounded int>0. 

Proof Let h eS satisfy h{0) = 1, h{x) = for x ^ B{0,e), an d h{x ) < 1 
for all X. Define r = min(t, TB(o,e))- Then, by Dynkin's formula ( |2.13[ ), 

Po{t > TBiO,e)) = 1 - Po{t < TBiO,s)) < h{0) - Eoh{Xr) 



Proof of Lemma \2.11\ Since there are few tools to handle functions non- 
vanishing at infinity, we use the definition of Ad- For e G (0,r], denote by 
D, the set of x such that dist(x,R'^ \ D) > e, and D'^ = R'^ \ D,. Our 
strategy is to show that away from the boundary of D (in Z^^) we have 
Pj^F P^f^ while near the boundary (in D \ D^), P^F ^ 0. 

Let Qeit) = Px{t > TB {x,e))- Notc that Qeit) docs uot depend on x G R'', 



and that by Proposition 2.12, ge{t)/t is a bounded function of t > 0. 



First, we consider x G -De, where e G (0,r] will be specified later. By 



Lemma 2.2, F G T>{A\ Cb). Since F vanishes on D'^, for x E we have 
\PtF{x)-P,''F{x)\ = \Ex{F{Xt)lt>rJ\ 

— \\F\\L°°{R.)Ex{lxteDAt>TD) ■ 

Observe that the condition Xf G and t > implies that before 
time t, the process first exits D, and after that it exits D'^. Furthermore, 
B{Xrj^,e) C D'j. a.s. Hence, by the strong Markov property, 

Ex{lxt£DAt>To) ^ Pxi^s first exists D and then B{Xrjj,e) before time t) 

= Ex{ge{t - TD)lt>Tn)- 
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Since Qe is increasing, and B{x,e) C D, we have 

E,{lx,eDAt>rn) < 9s{t)P.{t > Td) 

< 9s{t)P.{t > rBi,,s)) = {9eit))^- 

We conclude that for x E D^, 

PtF-F 



(2.15) 



< 



+ 



t 



|l°°(r) 



AF 



{9em 



L°°(R<*) 



t 



Now we consider x close to the boundary. Fix (small) rj > 0. By Corol- 
lary 2.6, AF is Lipschitz continuous. Hence, we can choose e G (0,r/2] 
small enough, so that |^F(?/)| < rj when dist(y, D \ D^) < e. Let (fn E S he 
the sequence approximating F as in the proof of Lemma 2.2 Let x E D\D^ 
and r = mm{t,To). By Dynkin's formula (2.13), 



E, 



Aipn{Xs)ds 



< E 



When t < TB{x,e)y we use the estimate 



\A(Pn{Xs)\ds <t sup \A(pn{y)\- 

y£B{x,e) 

When t > Tb{x,£)i we simply have 

\Aipn{Xs)\ds < i(:||^v'„||ioo(Rd). 

^0 

It follows that 

\E-,LPn{Xr) - ^n{x)\ < tP ^{t < TB(x,e)) SUp \A(fniy)\ 

y&B{x,e) 

+ tPx{t > TB{x,e))\\AiPn\\L^(Rd) 

<t\ sup \Afn{y)\ + ge{t)\\Afn\\L--iIl'') ■ 
\y(^B{x,e) J 

Recall that (fn and Aifn converge to F and AF, respectively, locally uni- 
formly, and that ||^(y9n||(2;2(Rd) < C||(/9„||f74(Rd) < C^||F||c'4(-j^d) for some 
C > 0. By taking the limit as n — t- oo and applying dominated conver- 
gence, we obtain that for a; G -D \ -D^, 



|E.F(X.)-F(x)|<t( sup \AFiy)\+C'gemF\\c.^^.^ 

y^B{x,e) 

We have |^F(?/)| < i] ioi y E B{x,e) (because dist{y,D \ D^) < e). Hence, 



E^F{Xr) - F{x) 



t 



AF{x] 



t 



<2r] + C^g, 
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Finally, -F(Xt-^) = a.s., and therefore 

E,F(X.) = E,{F{Xt)lt<rJ = PFF{^)- 

We conclude that 

PFF{x) - F{x) 



(2.16) 



t 



AF{x) 



< 27^ + Cg,{t)\\F\\ct(n^y 



Recall that as t — )■ 0"*", ge{t)/t is bounded. Furthermore, {PtF — F)/t 
converges uniformly to AF. Hence, (2.15) and (2.16) imply that for t small 
enough, 

P[^F - F 



t 



IdAF 



Since rj > was arbitrary, the proof is complete. 



□ 



2.4. Subordinate Brownian motions. From now on, we only consider 
one- dimensional Levy processes. The process Xt corresponding to /3 = 
1, 7 = and u vanishing, is the Brownian motion, running at twice the 
usual speed (that is, VarXj = 2t). We denote the transition density of the 
Brownian motion by ks{z), 



(2.17) 



ks{z) 



\/4 



TTS 



exp 



z 

'4^ 



s > 0, z G R. 



A Levy process Xt is said to be a subordinate Brownian motion if it can be 
written in the form Xt = Bz^, where Bg is the Brownian motion, Zt is a 
suhordinator (a non-decreasing Levy process), and Bg and Zt are indepen- 
dent processes. We assume that under P^;, Bs starts at x and Zt starts at 
a.s. To avoid trivialities, we assume that Xt is non-constant, that is, /3 > 
or V is non-zero. 

The subordinator Zt is completely described by its Laplace exponent ip{^): 
we have EqC-^^' = e-*^(«) (^ G C, Re^ > 0), where 



(2.18) 



Re^ > 0. 



Here /3 > is the same as in (2.10), and fi is the Levy measure oi Zf. a 



Radon measure on (0, cxd) satisfying m.m{l, s)fi{ds) < oo. The relation 
between the Levy measures n and u of Zt and Xt is 



(2.19) 



ks{z)fi{ds) 



zeR\{0} 



with ks{z) defined in (2.17). The Levy-Khintchine exponent of the subor- 
dinate Brownian motion Xt satisfies \I'(^) = ipi^"^)- 

The class of subordinate Brownian motions is naturally divided into com- 
pound Poisson processes (with bounded Laplace exponents ip) and processes 
corresponding to unbounded ip. The latter class consists of strong Feller pro- 
cesses: if Zt > a.s., then Xt = Bst has absolutely continuous distribution. 



This property is important for applications of Lemma 2.10 See |T0[ HH [66] 
for more information on subordinate Brownian motions. 
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Remark 2.13. One needs to distinguish the killed subordinate and the sub- 
ordinate killed processes. The first one is X/^ introduced above, while the 
other one is obtained by subordination of the killed Brownian motion. We 
emphasize that these two processes are essentially different. In particular, 
the spectral theory for the subordinate killed process is trivial: its eigenfunc- 
tions are the same as the eigenfunctions of the (insubordinate) killed Brow- 
nian motion, and only the corresponding eigenvalues are different. For a 
discussion of the relation between killed subordinate and subordinate killed 
processes, see e.g. j69]. □ 



We recall some standard definitions; see |66] for further information. A 
function / : (0, oo) — t- R is said to be completely monotone on (0, oo) if 
/(")(z) > for all z > and n = 0,1,2,.... A function / : (0, oo) ^ R 
is a Bernstein function if / is nonnegative and /' is completely monotone. 
By Bernstein's theorem, every completely monotone function is the Laplace 
transform of some Radon measure on [0, oo), and every Bernstein function 
has the form 

POO 

(2.20) f(z)=ci + C2Z+ {1- e-'')m{ds), z > 0, 

Jo 

for some Ci,C2 > and a Radon measure m on (0, oo) satisfying 
Jq mm{l, s)m{ds) < oo. Whenever we write /(O) when / is a Bernstein 
function, we mean lim2__>o+ /(^)- Note that if /(O) = 0, then / has the 



form (2.18); hence, a Bernstein function / is the Laplace exponent of a 



subordinator if and only if /(O) = 0. 

2.5. Complete Bernstein functions and Stieltjes functions. A Bern- 
stein function / for which the measure m in the representation ( |2.20[ ) has 
a completely monotone density function is said to be a complete Bernstein 
function (CBF in short; also known as operator monotone function). Re- 
call that Xt is a subordinate Brownian motion, /x is the Levy measure of 
the underlying subordinator, and ip is its Laplace exponent (so that the 
Levy-Khintchine exponent of Xt is iplC,"^)). 

Proposition 2.14. The following conditions are equivalent: 

(a) ip is a complete Bernstein function; 

(b) fi{ds) has a completely monotone density function fi{s); 

(c) for a Radon measure fiQ on (0, oo) such that the integral 

min(C"S C"^)Ato('^C) is finite, 

(2.21) m = /3^+' ^ 



(d) ip extends to a holomorphic function in C \ (— oo,0], which leaves 
the upper and the lower complex half-planes invariant; 

(e) v{z) is a completely monotone function on (0, oo). 

A symmetric Levy process satisfying (e) is automatically a subordinate 
Brownian motion. 
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Proof. Equivalence of (a) and (b) is just the definition of complete Bernstein 
functions. Equivalence of (b)-(d) is standard, see [66j, Chapters 6 and 11. 
Equivalence of (c) and (e) seems to be well-known to specialists (cf. ^62]), but 
a direct reference to the literature is difficult to find. The proof is a rather 
straightforward calculation involving Fourier transform and an application 
of Bernstein's representation theorem; for completeness, we provide the 
details. 

Suppose that (c) holds. Then, for ^ 

2, ,1 r i^oiio 



By an explicit calculation, for C > 0, 

2 2v^ 2 



cos 



Hence, by Fubini, 

= + ly - cos(e.)) f e-^l^l ^) dz. 

Comparing this with (2.10), we obtain 

(2.22) re-^W^^2(^, 

27r Jo VC 

which proves that ulz) is completely monotone on (0, oo). Conversely, if (e) 
holds, that is, if z/(z) is completely monotone on (0, oo), then, by Bernstein's 
theorem, i>{z) has the representation (2.22), and the above reasoning can be 
reversed to prove that Xt is a subordinate Brownian motion, and that (c) 
holds. □ 

The assumptions of the main theorems can be put in the following form. 

Assumption 2.15. The process Xt is a subordinate Brownian motion sat- 



isfying any of the equivalent conditions of Proposition 2.14 



We note the following immediate consequence of Proposition 2.14 

Corollary 2.16. The following conditions are equivalent: 

(a) f is a complete Bernstein function; that is, ( 2.20[ ) holds for some 
Ci,C2 > and a measure m with completely monotone density; 

(b) for some Ci,C2 > and a Radon measure mo on (0, oo) such that 

mm{s~^ , s~^)mo{ds) < oo, 

(2,23) /(^)=c,+c.. + i r^?^. .>0i 

TC Jo Z + S S 

(c) f{z) > for z > and either f is constant, or f extends to a 
holomorphic function in C \ (— oo, 0], which leaves the upper and the 
lower complex half-planes invariant; 
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(d) /(^^) — /(O) is the Levy-Khintchine exponent of a symmetric Levy 
process, whose Levy measure has a completely monotone density 
function on {0,00). □ 



The right-hand side of formula (2.23) defines a holomorphic function in 



C \ (—00, 0]. We often denote this extension by the same symbol /. We will 
also use the related notion of a Stieltjes function: / is said to be a Stieltjes 
function if there are Ci, C2 > and a Radon measure mo on (0, 00) satisfying 
Jq°° min(l, s~^)mo((is) < 00 such that 

c 1 1 

(2.24) f{z) = 5i + - + - / — rhoids), z > 0. 

z n Jq z + s 



Since l/{z + s) = e ^''e ^^dC,, we see that / has the representation (2.24) 
if and only if it is the Laplace transform of the measure 

Ci5o(rf0 + (£2 + 7i-^Cmo{0)d^ 

on [0,00), where 60 is the Dirac delta measure. In other words, / is the 
Laplace transform of a completely monotone function, plus a nonnegative 
constant. 

We list some standard properties of complete Bernstein and Stieltjes func- 
tions. 

Proposition 2.17 (Chapter 3 in [37] and Chapter 7 in [SB]). Suppose that 
f and g are not constantly equal to 0. 

(a) The following conditions are equivalent: f{z) is a CBF; z/f{z) is a 
CBF; l/f{z) is a Stieltjes function; f{z)/z is a Stieltjes function. 

(b) If f, g are CBF, a,C > and a e (0, 1), then f{z) + g{z), Cf{z), 
f{Cz), {z — a)/{f{z) — /(a)) (extended continuously at z = a), 
f lg{z)) and (/(2"))^/" are CBF. □ 

One of the fundamental properties of complete Bernstein and Stieltjes 



functions is that the characteristics in representations (2.23) and (2.24) 
can be easily recovered from (the holomorphic extension of) the function. 
We illustrate this for Stieltjes functions. By 5q we denote the Dirac delta 
measure at 0. 

Proposition 2.18. Let f he (the holomorphic extension of) a Stieltjes func- 
tion with representation ( |2.24[ ). Then 

(2.25) ci = lim f{z), C2 = lim {zf{z)), 

z-s>oo 2-S>0+ 

and 

(2.26) niQids) + TTC2So(ds) = lim (—Im f(—s + ie)ds), 

£-S>0+ 

with the limit understood in the sense of weak convergence of measures. 



Formula (2.26) is well-known in complex analysis, it is a variant of the 
Sokhotskyi-Plemelj formula for the Cauchy-Stieltjes transform in the upper 
complex half-plane. The proof of a version for complete Bernstein func- 
tions is contained in Theorem 6.2 in |l66j (see Corollary 6.3 and Remark 6.4 
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therein). We provide a simple argument based on the representation theo- 
rem for harmonic functions in half-plane. 



Proof of Proposition 2.18 Since — Im / is a nonnegative harmonic function 
in the upper complex half-plane, it can be represented using the Poisson 
kernel, 

r>00 



(2.27) ~lmf{x + iy) 

for unique costant C = lim 



Cy 



y 



-oo y^ + {x- s) 



m{ds) 



— Im f{iy) I y) and unique measure m[ds) 



lim£_j.o+(— Im/(s + ie)ds). On the other hand, by (2.24), 



Im /(x -|- iy) 



C2y 



y^ + 



y^ + [x + s) 



mo{ds). 



Formula (2.26) follows simply by comparing the above two representations. 
The expressions (2.25) are direct consequences of (2.24) and dominated 
convergence. □ 

We need three more properties of complete Bernstein functions. 

Proposition 2.19. If f is a complete Bernstein function and a > 0, then 
g{z) = {l — z/a)/{l — f{z)/f{a)) and h{z) = logg{z) are complete Bernstein 
functions. 

Here we extend g and h continuously ai z = a, g{a) = f{a)/{af'{a)) and 
h{a) = log(/(a)/(af (a))). 

Proof. By Proposition 2.17K b), g(z) = (f{a)/a)(x — a)/{f(z) — f{a)) is 
a CBF. Note that (7(0) = 1, and so also g{z) — 1 is a CBF. Furthermore, 
log(l-Fz) is a CBF. Again by Proposition |2.17K b), h{z) = log{l + {g{z) -1)) 
is a CBF. □ 

Lemma 2.20. Let f be (the holomorphic extension of) a complete Bernstein 
function, and let a > 0. Define 



1 f fjia) _ fj-ia) 
2ia \z — ia z + ia 
2 I ^2 



z E C \ {—ia, ia}. 



Then g{z) = f*{z) — f{z)/(z + a ) is a Stieltjes function, and it is the 
Laplace transform of a completely monotone function G on (0, 00). Fur- 
thermore, G is the Laplace transform of a finite measure on (0, 00). 



Proof. Let / have the representation (2.23), and let h{z) = f{z) — Ci — C2Z. 
Fix z e C \ (—00, 0]. Since 



z + s a^ + 



la 



2ia{z 



ia) s + ia 
1 



-la 



we have 



mo{ds) 



ia){z + ia) s 
h{ia) 



la 



a^ + z^ z + s 



z + s a"^ + s"^ 



1 
2ia 



la 



h{—ia) 
z + ia 



h{z) 



+ z"^ 
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Furthermore, 

1 Ci + C2ia Ci — C2ia\ C\ + C2Z ^ 



lia \ z — ia z + ia J a'^ + z'^ 
We conclude that 

1 p 1 moids) 1 //(za) /(-za)\ f{z) 



TT Jo 2; + s + 2ia \z — ia z + ia J a^ + z^ 
By the definition, g is the Stieltjes function, and by the remark follow- 
ing (2.24), g is the Laplace transform of a completely monotone function 
G{x) = C'y{x), where •y^ds) = (l/7r)(a^ + s^)~^mo((is) is a finite measure 
on (0, 00). □ 

Proposition 2.21. /// is a complete Bernstein function, then 

(a) < zf'{z) < f{z) for z > 0; 

(b) < -zf"{z) < 2f'{z) for z > 0; 

(c) \f{z)\ < (sin(£/2))-V(kl) forze C, \ ATgz\ < n - e, e e (0,7r); 

(d) \zf'{z)\ < (sin(£/2))-i/(kl) for z as m (c); 

(e) \f{z)\ < C{f,e){l + \z\) for z as m (c). 

In (c)-(e), the holomorphic extension of f is denoted by the same symbol. 



Proof. All statements follow easily from the representation (2.23) and dif- 
ferentiation under the integral sign. Indeed, (a) reduces to 

1 sz mo(ds) 1 z mo(ds) 

0<C2Z+- . , <C, + C2Z+- — , 

tvJq {s + z)^ s tvJq s + z s 

which follows from the inequality sz/{s + zY < ^/{^ ^ ^) fo^' s^z > Q. In 
a similar manner, (b) follows from 2sz/{s + zY < 2s/{s + z^. To prove 
(c)-(e), observe that for 2; G C with Argz G [—it + e,TT — e], 

\s + z\^ = + 2sRez+ \z\^ > - 2s\z\ cos{e) + 

Hence, by a simple calculation, 

\s + z\'> i^^^M (, + 1^1)2 = (sin(e/2))2(. + 1^1)1 

Therefore, (c) is a consequence of \z/{s + z)\ < (sin(£:/2))~-'- |2;|/(s + \z\), 
and (d) is proved using \sz/{s + zY\ < \z\/\s + z\ < (sin(£:/2))^-^|2;|/(s+ \z\). 
Finally, to prove (e), use (c) and the inequality < (l + |z|)/(s+l), 

verified easily by a direct calculation. □ 

For further properties of complete Bernstein functions and related no- 
tions, see Chapters 6 and 7 in [66j, or Chapter 3 in |37) . 

3. Transformation of the Laplace exponent 

Recall that Xt = is a subordinate Brownian motion, ipiO the 
Laplace exponent of the subordinator Zt, and \E'(^) = is the Levy- 

Khintchine exponent of Xj. Assumption 2.15 is in force; that is, ip is as- 
sumed to be a complete Bernstein function. In this section we introduce 
the operation ip tp"^ . Although it will be used for complete Bernstein 
functions only, the definition is quite general. 
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Definition 3.1. When ^ is a positive function on (0, oo), and tlie integral 
Jq°° min(l, ^^^) \ogip{^)d^ is finite, then we define 

(3.1) = exp I - J ^ ) ' ^ > 0- 

Proposition 3.2. For -ip as Definition \3. 1\ 



(3.2) *'tt)=exp|ijJ"iMidC), 5>0. 



Proof. Substitute ,^(^ = s in (3.1 ). □ 



Note that (3.2) defines a holomorphic extension of ip^ to the half-plane 



Re^ > 0. We denote this extension by the same symbol ip'^iC,)- 

Remark 3.3. The definition of ip"^ lies at the very heart of the Wiener- 
Hopf method. Roughly speaking, if ^/'^(^) can be extended continuously to 
the boundary of the region Re^ > 0, that is, to the imaginary axis, then 
we have = ilj\i^)'ip\—i^) , a Wiener-Hopf factorization of For 

complete Bernstein functions, this is formally proved in Lemma [3.8| below. 
The Wiener-Hopf method is described in Section |4j 

The function ip"^ plays an important role in fluctuation theory of Levy 
processes. When Xf is a symmetric Levy process, then its Levy-Khintchine 
exponent \E' is symmetric, ^'(0 = V'(^^) fo^' some nonnegative ip on [0,oo). 
Suppose that Xt is not a compound Poisson process. By Corollary 9.7 
in [28j, the Laplace exponent k,{z,C,) of the bivariate ascending ladder pro- 
cess satisfies 

^{z,O = {z + 0{O, z>0,^>0. 
This relation was used extensively in |l3l HH [55]. □ 



The fundamental Lemma 3.8 states that if -0 is a CBF, then is a CBF. 



Before we prove it, we establish some simple properties of ip"^ . 

Proposition 3.4. Whenever both sides of the following identities make 
sense, we have: 

(a) {1/0 = (7/'")t = (^t)a ^ {i)xi)2)'^ = ^M; 

(b) {C^tpy = CtP^ for C > 0; 

(c) when pj{^) = ^, then ^p\^) = ^; 

(d) if appropriate limits of ip exist, then the corresponding limits of 
ip'^ exist, and limg^o+^^(0 = (lim^^o+ ^(0)^^^ limg^oo^"^(0 = 



Proof. Properties (a) follow immediately from formula (3.2). In (b), one 



uses also /^^^/(^^ + C)dC = vr/2. For ^^(0 = ^, by ^ we have 

1 /■°° 21oge + 2 1ogC 



^t(e) = exp(-/ 



exp(ioge + - r^^Cl =e; 



SUBORDINATE BROWNIAN MOTIONS IN HALF-LINE 25 

here we used the identity /q°°(1 + s"^)"^ log s ds = 0, proved easily by a 
substitution r = 1/s. Finally, to show (d), assume that a finite, posi- 
tive limit ip{0) = limg_>o+ "^(0 exists. It can be proved that the functions 
log'i/'('^^C^)/(l + C^) of C ^I'e uniformly integrable as ,^ — )■ 0+ (this is clear 
when ip is Si CBF; we omit the details in the general case). It follows that 

riogv-m /-logM .iogv<o) 



and hence limg_^o+ V^^(0 = a/V'(0)- When the limit of is or oo, one uses 
in a similar way Fatou's lemma; we omit the details. Limits at oo are dealt 
with in a similar way. □ 



Proposition 3.5. Let ip he as in Definition \3.1\ If Ci,C2 > 0, a G R 
and ip{C,) < {Cf + C|^)" for all > 0, then, with the same constants, 
< |Ci + <^2^l"; Re^ > 0. In a similar manner, ifipi^ > (C^ + Cl^)" 
for^ >0 for some Ci,C2 > 0, a G R, then |^t(^)| > ic^ + CsCl", Re^ > 0. 

Proof. For ^ = t + is (t > 0, s G R), we have 
l^^'m' = exp r Re log<;'(C')<ic) 

= '''' (; r + (C - s)^ ^ t^ + (C + .)0 '"'"''^'''^ 

Suppose that ^/'(C) < (Cf + C|C)" for C > 0. Then 

log^(C') < alog(C2 + Cie) = 2aRe(log(Ci - C^tQ), 

and the right-hand side extends to a bounded below harmonic function in 
the upper half-plane Im^ > (with vanishing linear term in the represen- 
tation similar to (2.27)). By the Poisson formula, 

I f°° t 

Re(log(Ci - C2<))rfC = Re(log(Ci - C2i{s + it))). 



TT t^ + iC- S) 

Hence, 

|^t(^)|2 < exp {2a Re(log(Ci - C2i{s + it)))) 

= ((Ci + Cat)' + Cls^)"" = |Ci + Cs^r" 

This proves the first statement. The other one follows by reversing the 
inequalities in the above argument. □ 

Remark 3.6. In Proposition 2.1 in [55] it is proved that if ip{^) and ^/^(O 
are increasing, then iIj{C,)/2 < ■?/'^(^) < 2'iIj{^). For complete Bernstein 
functions, this was independently proved in Proposition 3.7 in [44]. □ 
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Recall that under Assumption 2.15 ip is a. complete Bernstein function. 
Hence, ip extends to a holomorphic function in C\(— oo, 0]; we use the same 
symbol ip for this extension. First, we define an auxiliary function 

1 riog^pi^e) 



(3.3) fiO = exp 



TT 



eGC\(-oo,0], 



so that = /(^^). Note that (/(O)^ is an (integral-type) weighted 

geometric mean of the family of complete Bernstein functions V^(^C^) (C > 
0), so that (/(O)^ is a complete Bernstein function. This is formally proved 
in the following simple result. 

Lemma 3.7. If ip is a complete Bernstein function, then f{C,) and (/(O)^ 
are complete Bernstein functions. 



Proof. Clearly, /(^) > for ^ > 0. Furthermore, 
1 /"^ Imlog^(ea 1 

When Im^ > 0, we have ATgip{C,('^) 



Arg/(0 



dC 



TT 



Arg^(eC' 



dC. 



1 + 
(0,7r). It follows that Arg/(0 G 
(0,7r/2), that is, Im/(0 > and Im(/(0)' > 0. Finally, fiO = M, 
so Im/(^) < and Im(/(^))^ < when Im^ < 0. The result follows by 
Corollary [2l6| ^c). □ 

Below we prove that ip'^lC,) = /(■C^) is a complete Bernstein function; this 
is stronger than the assertion of Lemma 3.7, see Proposition 2.17[ b). 

Lemma 3.8. If is a complete Bernstein function, then ip'^ is a complete 
Bernstein function. The holomorphic extension of ip'^ ( denoted by the same 
symbol) satisfies 



(3.4) 



eeC\R. 



Remark 3.9. The first statement of the above lemma was proved inde- 
pendently in Proposition 2.4 in ^3] using Theorem 6.10 in [66j. After the 
preliminary version of the article was made available, Jacek Malecki pointed 



out to the author that the other statement of the lemma, formula (3.4), at 



least for purely imaginary ^, can be deduced from the Wiener-Hopf fac- 
torization of the Levy-Khintchine exponent in fluctuation theory, see for- 
mula (VI. 4) in [5j. The novelty of Lemma 3.8 lies in the combination of 



the two parts, which is one of the key steps in the derivation of the explicit 
formula for the eigenfunctions Fx in Theorem |1.1[ □ 



Proof. Let / be (the holomorphic extension of) the complete Bernstein func 



tion defined in (3.3) (see Lemma 3.7); hence, V^(0 — /(^^) fo'^ ^ > 



and this identity defines a holomorphic extension of ip"^ to the half-plane 



Re^ > 0. When Re^ > 0, by (|3.3|) we have 

1 



^Ho = fie 



exp 



TT 



dC 



10 i+c^ 

For a fixed ^, the integrand on the right-hand side is a meromorphic function 
of ( in the region Re(^C) > with a simple pole at —i. Furthermore, 
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(a) 



\ 
(b) 



Figure 1. Two contours used in the proof of Lemma 3.8 



Dashed hne depicts the boundary of the region in which the 
integrand is meromorphic. In (a), Re^ > and Im^ > 0, 
and the integrand logip{^'^('^)/{l + C^) is meromorphic in the 
region Re(^C) > with pole at —i, outside the contour of 
integration. In (b), Re^ < and Im^ > 0, and the integrand 
log'?/'(^^C^)/(l + (^) is meromorphic in Re(— ^C) > with a 
simple pole at i. 



2.21[e), for any e > the integrand decays at least as 



by Proposition 

fast as log ICI^ when |C| — )■ oo in the region —tc/2 + e < Arg(^C) < 
7t/2 — 6. Hence, by a standard contour integration (using contours shown in 
Figure [T] (a) with r — )■ oo; the pole at —i is outside the contour) and then 
parametrization C = s/y^ of [0, ^"^''^oo), we obtain that 

-V2oo 

^t(e)=exp|3 / ^^v:\: ' dc 



exp 



TT 



TT 



ds 



The formula on the right-hand side clearly defines a holomorphic function 
in C \ (—00,0] (integrability follows by Proposition 2.2lK e)), which is the 
holomorphic extension of V^^(0- Suppose that Re^ < and Im^ > 0. Using 



the parametrization C, 

= 




""^/^oo), we obtain that 



dc 



For a fixed ^, the integrand on the right-hand side is a meromorphic function 
of ( in the region Re(— ^C) > with a simple pole at i. Furthermore, again 
by Proposition 2.21[ e), for any e > 0, the integrand decays at least as fast 
as C~^log|CP when \(\ — )■ oo in the region —n/2 + e < Aig{—^Q < 7r/2 — e. 
Hence, again using standard contour integration (along contours shown in 
Figure [l] (b) with r — )■ oo; this time the pole at i is inside the contour) and 
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the residue theorem, we obtain that 

1 riog^iee) 



^"^(0 = exp 



exp 



TT 



TT 







d( + 2i Resj 



2/-2\ 



fie 



Recall that Re^ < and Im^ > 0. Hence Re(-^) > 0, so that /(^^ 
/((-O^) = i^H-O- Formula ([3^ follows when Re^ < and Im^ > 0. 



Since ip and / are complete Bernstein functions, by Corollary 2.16K c), 
hi{^) = Arg ■?/;(— and /i2(0 — A^S/(^^) ^ire bounded harmonic func- 
tions in the region Re^ < 0, Im^ > 0, taking values in [0,7r] and [— 7r,0], 
respectively. Hence, hi and /12 are Poisson integrals of their (bounded, mea- 
surable) boundary values gi, g2 (defined on (—00, 0]U[0, too)). Furthermore, 
we have gi{is) = Arg-?/;(s^) = and g2{—s) = Arg/(s^) = for s > 0, and 
gi{—s) G [0,7r], g2{is) G [— 7r,0] for almost all s > 0. 
Observe that when Re^ < 0, Im^ > 0, we have 

Arg^t(^) = ArgV^(-e') - Arg/(e') = hiO ~ /^2(0- 

Hence, /i(^) = Arg (^) is a bounded harmonic function in the region 
Re^ < 0, Im^ > 0, with boundary value g given hj g = gi on (— oo,0], 
and g = — (?2 on [0,200). It follows that g{C,) G [0, vr] for almost every 
boundary point ^, and so h{^) G [0, vr] for all ^ such that Re^ < 0, Im^ > 0. 
It follows that Im^/'^(^) > in this region. We have already seen that 
Imip^^) = Im/(^2) > when Re^ > 0, Im^ > 0. 
and so Imip'^l^) < when Im,^ < 0. By Corollary |2.16 



and (3.4) extends to all ^ G C \ R. 



Finally, tfj^iO = ^KO, 
T/^t is a CBF, 
□ 



We note a third formula for (the holomorphic extension of) ^^(0, other 
than (3.1) and (3.2): when ip is a complete Bernstein function, then 



(3.5) ,.,0^exp(if 



^gC\(-oo,0]; 



see the above proof of Lemma 3.8 

In Section |4| the transformation ip ip'^ is applied to the function 



(3.6) 



MO 



1 - e/A^ 



A > 0, ^ G C \ (-00,0). 



i-V'(0MA2)' 

This definition is continuously extended at ^ = A^ by letting ipxlX'^) 
^(A2)/(AV(A2)). We denote ^1(0 = (^a)^^- 



By Proposition 2.19, ipx and logV'A are complete Bernstein functions, and 



iJx{0) = 1. It is easy to see that ^^(A) = {i:{X^)\ij"{\^)\)/{2\'^{f{y))^) 
(note that ■?/;"( A^) < 0). Since ip is concave, its graph lies below the tangent 
line at ^ = A^, that is. 



(3.7) 



MO < 



AV(A2) 2A2(^'(A2))2 



A,e>0. 
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Formula (3.7) (plus ^ — < ^) combined with Proposition 3.5 yields that 




^(A2)|^"(A2)|^ 



A > 0, Re^ > 0. 



(3.8) ^1(0 < ■ Y 2\\i,'{\^)Y 

We remark that more detailed estimates of ipx and ip^ can be found in [56]. 

4. Derivation of the formula for eigenfunctions 
Below we prove Theorem Our strategy is as follows. First we state a 



distributional version of the spectral problem for ^(o,oo), and in Lemma 4.2 
we prove that a distributional solution is automatically a strong solution. 
The (simpler) case of compound Poisson processes requires a different ap- 



proach and is studied separately in Proposition |4.4[ Next, we rephrase the 
distributional problem as a Wiener-Hopf equation, and we solve it using 
the Wiener-Hopf method. Finally, we combine the two parts and prove 
Theorem In addition, we establish some basic properties of the eigen- 
functions. 

4.1. Distributional and strong eigenfunctions. In this part, Xt can be 
an arbitrary one-dimensional symmetric Levy process with Levy-Khintchine 
exponent \E'. We define by the formula '^{C) = '^iC'^)- Note that 
need not be a Bernstein function. 

Recall that A is the generator of Xf, and ^(o,oo) is the generator of the 
process Xf killed upon leaving (0,oo), defined on an appropriate domain. 
Here we only consider the Cb and Cq generators, formal definitions 
are given in Preliminaries. Recall also that A E S' is the distributional 



generator of Xt, see Definition [2^ The Fourier transform of the distribution 
A is —ipi^"^), and since A is symmetric, it is equal to its reflection A. 

Let F G L°°((0, oo)). As usual, F is extended to entire real line by setting 
F{x) = for X < 0. Given some regularity of Xf, if F G Co((0, oo)) then. 



by Lemma 2.10, the condition ^(0,00)-^ = —Tp{X )F is equivalent to the 
apparently weaker condition A * F{x) = —tp{X'^)F{x) for x G (0, 00). This 
motivates the following definition. 

Definition 4.1. Let A > 0. A tempered distribution F is said to be a dis- 
tributional eigenfunction of ^(0.00), corresponding to the eigenvalue —iplX'^), 
if F is supported in [0, 00), iS'-convolvable with A, and A* F + ip{\^)F is 
supported in (— oo,0]. Informally, 

F = on (-00,0), 

f4 1) 

^ ' ' A*F = i){X^)F on (0,00). 

Lemma 4.2. Suppose that Xt is a symmetric Levy process, which has the 
strong Feller property, and that is a regular boundary point of (0, 00). 
Let A > 0, and suppose that F is a distributional eigenfunction of A(o^oo), 
corresponding to the eigenvalue —'ip{X'^). If F is bounded and continuous on 
R, and if 

lim (F(x) - C sin(Ax + ^)) = 
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for some C,^ eR, then F G ^'(^(o.oo); ^°°) and A{o,oo)F = -ip{X^)F. 



Proof. Since F does not vanish at infinity, we cannot apply Lemma 2.10 



directly to F to prove that F G ^^(^(0,00)! L'^)- Thus, we decompose F into 
the sum of a Co((0, 00)) function and a smooth bounded function. 



Denote F*{x) = Csin(Ax + ^). By Remark 2.1, F* G V{A; 1"=^) and 



AF* = -iIj{\^)F*. By Proposition [Z5| A* F* = -ip{X^)F* . We choose an 
infinitely smooth function h such that h{x) = for x < 1 and h{x) = 1 for 
X large enough. Furthermore, we require that 

/oo 
h{z)F*{z)v{z)dz = 0. 

(Note that such a choice is always possible, as we do not assume that < 
h{x) < 1.) We write F = fi + /a, where /i = hF* and f^ = F - hF*. 
The function /i is infinitely smooth, with all derivatives bounded, so by 



Lemma 2.2, /i G V{A]Cb). We have fi{x) = for a; < 1, and by (2.8) 



and (4.2), ^/i(0) = 0. Hence, by Lemma |2.11| and Proposition 2.5, we 
have fi G 'D{A(o,oo); L°°) and ^(o,oo)/i(a;) = Afi{x) = A* fi{x) for x G 
(0, 00). We claim that also /2 G T>{A(^o^oo)] L°°) and A(^o,oo)f2{x) = ^4 * f2{x) 
for X G (0,00). Once this is proved, we have F G 1^(^(0,00);-^°°) and 
A(o^oc)F'{x) = A * F{x) = —ip{X^)F{x) for a; > 0, as desired. Hence, it 
remains to prove the claim. 

We will apply Lemma 2.10 to f2 = F — hF*. Since F and h vanish on 
(-00,0], we have /2(x) = for x < 0. Also, f^ = {F - F*) + (1 - h)F*, 
and both F — F*, {1 — h) vanish at infinity, so that lim^.„>oo f2{x) = 0. In 
other words, /2 G Co((0, 00)). Furthermore, 

A*f2 = A*F-A* (hF*) = A* F - A* fi. 

By assumption, A*F restricted to (0, 00) is equal to —ip{X'^)F{x). We have 
already seen that A*fi = Afi is a continuous function on R which vanishes 
at 0. Hence, A * f2 restricted to (0, 00) is equal to a continuous function, 
vanishing continuously at 0. Furthermore, for x > 0, 

A * /2(x) =A* F{x) -A* F*{x) + A*{{1- h)F*){x) 

= -^{X^){F{x) - F*{x)) + A*{{1- h)F*){x). 

By assumption, the first term on the right-hand side converges to as x — )■ 
00. Since 1 — h{x) = for x large enough, also A* {{1 — h)F*){x) vanishes 
as X — )■ cxD (by Lemma 2.2 and Proposition |2.5 ). We conclude that A* f2 
restricted to (0, 00) is equal to a Co((0, 00)) function. By Lemma 2.10, /2 G 
^^(Ao.oo); C'o) (and therefore also /2 G V{A(o,ooy, L°°)), and ^(o,oo)/2(a;) = 
A * /2(x) for X > 0. □ 

Remark 4.3. Continuity of F at is essential. Indeed, when Xt is the 
Brownian motion (i.e. ip{C,) = ^), then F(x) = cos(Ax)l(o,oo)(a^) is the 
distributional eigenfunction of ^ = d'^/dx^ in (0, 00), but F is not in the 
L°°(R'^) domain of ^(0,00) • In particular pI^^'°°^F is not equal to e~'^^*F, as 
can be verified by a direct calculation. □ 

When Xt is a compound Poisson process, the proof is much easier. 
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Proposition 4.4. Suppose that Xt is a symmetric compound Poisson pro- 
cess. Let A > 0, and suppose that F G L°^((0, oo)) is a distributional 
eigenfunction of ^(o,oo); corresponding to the eigenvalue —iplX"^). Then 
Ao,oo)F = -V^(A2)F.' 



Proof. The result is a straightforward application of Lemma 2.7 and the 



fact that A is the signed measure ^{dz) — u{lV)6Q{dz). □ 
4.2. Wiener-Hopf factorization. A method of solving problems of the 



form (4.1) was found in [72j in the case when A is an integrable function 
satisfying some growth conditions. For this reason, problems of this kind are 
called Wiener-Hopf equations, and the algorithm for finding their solutions 
is the Wiener-Hopf method; see [3T| [i6] for a detailed exposition. We begin 
with a brief description of the method. A version of Paley- Wiener theorem 
states that a function is supported in [0, oo) or (— cxd, 0] if and only if its 
Fourier transform extends holomorphically to the upper or lower complex 



half-plane, respectively. Therefore, (4.1) can be re- written as 



J^F{C,) extends to the upper complex half-plane, 

(C — ■ip{^'^))J-'F(^) extends to the lower complex half-plane, 

where C = ipiX"^)- The main ingredient of the Wiener-Hopf method is 
factorization of the symbol C — into the product of two functions, 

one extending holomorphically to the upper half-plane, the other one to the 
lower half-plane (the Wiener-Hopf factors). Then J^F is taken to be the 
multiplicative inverse of the first factor. 

The fundamental condition for the Wiener-Hopf method is that the sym- 
bol does not vanish. In our case, we have C — ip{^^) = for ^ = ±A. For 
this reason, we apply the Wiener-Hopf method to the regularized symbol 

i-^(a/v^(A^) _ 1 
i-e/A2 Me)' 



here we use the notation of (3.6). 



From now on. Assumption 2.15| is in force. It is required in order to 
establish some properties of 1/'?/'a(C^) and its Wiener-Hopf factors. In par- 
ticular, since ip\ is a complete Bernstein function, the factorization is given 



by Lemma 3.8 



The Wiener-Hopf factors are ip\{—iO aiid ip^^iC,), and the solution is ex- 
pected to satisfy J^F{^) = V^5(-iO/(l - ^VA^), up to a multiplicative 
constant. Due to the singularity in the denominator, we prefer to work 
with the Laplace transform of F. Below we provide details for the above 
idea. 

Fix A > 0, and define a tempered distribution 

A^ = ^P{X')5o + A, 
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where A is the distributional generator of Xt and 60 is the Dirac delta 
measure. Hence, J^A{^) = ipiX"^) — We decompose Ax into the 

convolution of the distribution 

(4.3) 5 = ^ {X% + S'^) 

(the singular part) and two infinitely divisible measures i?+ and i?_ {regular 
parts) supported in [0, 00) and (—00, 0] respectively. 
Clearly, 



(4.4) :FSiO = ^ (A^ - e) = ^(A^) (1 - I 



We define the tempered distribution R by means of Fourier transform, 

with ipx defined in ( |3.6[ ). Since S is compactly supported, S and R are 
5'-convolvable, and since J^S{^)J='R{^) = ^(A^) - ^(^^^ = J'Ax, by the 
exchange formula we have 

Ax = S*R. 

Lemma 4.5. There is a probability measure R+ supported in [0,oo), such 
that (see ) 

Furthermore, i?+ /ias a completely monotone density on (0, 00) (but it may 
have an atom at 0). 

If R-{E) = Rj^{—E) is the reflection of then R = i?+ * i?_ and 



Proo f. By Lemma 3.8, ip^ is a complete Bernstein function. By Proposi- 
tion 2.17 a), l/ipl is a Stieltjes function, and therefore it is the Laplace 
transform of a measure whi ch h as completely monotone density func- 
tion on (0,00). By Proposition Qd), 1/V^1(0) = l/ipx{0) = 1. Hence, 
i?+ is a probability measure. Fina 

holomorphic extension of ip\ is denoted by the same symbol). 
Let i?_ be the reflection of R^. By (3.8), 



ly, J^R+iO = lM(-^0 (where the 



J^R(^) = = = J^R_^{^)J^R_{^) 

Me) i^iHOi^lm 

for ,^ G R. The lemma follows by the exchange formula. □ 

Lemma 4.6. Suppose that F & S' is a distribution supported in [0, 00), for 
which F * if is bounded for any ip & S. If F * {S * R+) is supported in {0}, 
then F is a distributional eigenfunction of ^(0,00) ■ 
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Proof. Note that both F and S * R+ are supported in [0, oo), and so they 
are iS'-convolvable. Furthermore, for any (pi,(p2,'P3 G 5, F *(pi is bounded, 
and {S*Rj^)*!f2, R-*<f3 are integrable. Therefore, the iS'-convolution of F, 
S*R^ and R_ is associative. In other words, if we let F = F*{S*R^), then 
F*R^ = Ax*F. Clearly, if the support of F is {0}, then Ax*F = F*R_ 
is supported in (— oo,0]. □ 

4.3. Formula for eigenfunctions. Lemma |4.6| enables us to describe a 



family of solutions to (4.1). If F * {S * R^) is supported in {0}, then 
J^F-J-'{S*R^) = P for some polynomial P. Therefore, the Fourier transform 
of a distributional eigenfunction F is expected to have the form P/J^{S * 
R+); however, extra care is needed because of the zeros of the denominator 
J^(^ * R+) at ±A. 



Let (see ( |46| ) 



(4.7) f(n = = ^^^tMl ^ > 0, 

where 

(4.8) cx = 



'AV(A2 



^a(A2) V ^(A2 

is a normalization constant. Since is a complete Bernstein function, / 
extends to a meromorphic function in C \ (— oo, 0], with poles at ±iX. We 



denote this extension again by /. By (3.2) and Proposition |3.4K b) 
(4.9) 



^1(0 = VM>^ exp ( - r log dC 



'^''^ exprif^logW^,, 



" v^io e + e " ^(A2)-^(c^ 

this holds true for complex ^ with Re^ > 0. It follows that f{C,) is equal 
to the right-hand side of (1.1). In order to prove Theorem |1.1| we need to 



show that / is the Laplace transform of a function F of the form ( |1.2[ ), and 
that F is a (distributional, and hence L°°) eigenfunction of ^(o,oo)- 
We define 

(4.10) r(0 = ^($^-4^]. (eC\l-^X,^X} 



and 

(4.11) gio = no-fio 



2iX \ ^-iX ^ + iX 

1 fcx^itX) cx4{-tX)\ cx^iO 



2tX\ ^-iX ^ + iX / A2 + ^2 



Note that Ci(^ — 26*2) ^ is the Laplace transform of Cie *'"^^l(o,oo)(2;) 
Hence, /* is the Laplace transform of 

(4.12) f-W^*'^i^:^^fIt^^i,„,.,W. xeR. 
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By Lemma 



2.20 



(applied to a CBF c\ip\), g is a. Stieltjes function, it is the 
Laplace transform of a completely monotone function G{x), and G is the 
Laplace transform of a finite measure. Hence, / is the Laplace transform of 
F{x) = F*{x)-G{x). 

Proposition 4.7. We have ip^^iiX) 
by g. 

Proof. By ( [49| ), for e G (0, A) we have 

1 



2\ g±ii9 



^, where -t^x is given 



exp 



TT 



iX + e 



log , 



We claim that the integrand on the right-hand side is dominated by an 
integrable function. Then, by dominated convergence. 



lim ipxi'^X + e) 

£-5-0+ 



exp 



iX 



1 ^^(^')^^ 



as desired (see (1.3) and (4.9)). Hence, it remains to prove the claim. 
Observe that 



iX + e 






iX + e 




{tx + ey + e 


t{x-0- 




\i{X + C)+e\ 



< 



X + e 



< 



2X 



|A-C|(A + C) - |A-C|(A + C)' 
Since ipx is nonnegative and concave, we have i'\{C'^)/ip\{X^) > C^/X^ for 
C < A, and '?/'a(C^)/V^a(A^) < C'^/X^ for C > A. This and monotonicity of ipx 
gives that 



log 



It follows that 



^a(A2) 
iA + e 

(2A + £)2 + C2 

which is integrable in C > 0. 



< 



log 



log 



^a(A2 



< 



2|logC-logA|. 

4A| logC - logA| 
IC-A|(C + A) 



□ 



The above proposition and formulae (4.8) and (4.12) give that F*[x) = 
sin(Aa; + d\)l{Q^ao){x)- We summarize our development in the following 
result. 

Corollary 4.8. The function f is the Laplace transform of a function 
F{x) = F*{x) — G{x), where 

F*{x) = sm{Xx + ?9A)l(o,oo)(a;), 

G is completely monotone on (0, oo) and it is the Laplace transform of a 
finite measure on (0,oo). □ 

Lemma 4.9. The function F constructed above is a distributional eigen- 
function of A(o,oo)- 
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Proof. By Q, we have CS = {ij{X^)/X'^){X'^ + ^^). Since F = F* - G 
and both F* and G are bounded, also F is a bounded function. Hence, 
F is 5'-convolvable with S * R+, and ■ {CS{OCR+{0) = fiO ■ 

((^(A2)/A2)(A2 + ^2)/^t^^)) _ ca^(A2)/A2 when Re^ > 0. By the exchange 
formula, F * {S * = {cxipi^X"^)/ X'^)So, which is a distribution supported 
in {0}. An application of Lemma 4.6 completes the proof. 



□ 



Remark 4.10. Clearly, also the distributional derivatives (of all orders) of 
F are distributional eigenfunctions of .4.(o,oo)- However, they are not in the 
L°^((0, oo)) domain of ^(o,oo)- D 



In order to apply Lemma |4.2[ we need to show that F vanishes continu- 
ously at 0. 

Lemma 4.11. The function F has a nonnegative right limit at 0. If is 
unbounded, then F vanishes continuously at 0. 

Proof. Since G is non-increasing (even completely monotone), the limit of 
F{x) as X —7- 0^ exists. Since CF = /, where / is given by (4.7), we have 



lim F(x) 

a; -5-0+ 



hm ^/ = lim 

5-s>oo 5^oo -|- 




A2 + e e 



The above expression is nonnegative, proving the first part 



3.4 a,c) 



Furthermore, 
If ip is unbounded on [0, oo) 



by Proposition^ 

then '0a(O/^ tends to as ^ — )• oo (see ( |3.6 )). Hence, by Proposition 3.4[ d), 
i^xiO/^ converges to as ,^ — )■ oo, and therefore lima,_^o+ P{x) = 0. □ 



Proof of Theorem \l . 1 In the remark following (4.9) it was already observed 
that CF{C,) = /(O coincides with the right-hand side of the formula (1.1) 
for CFx, so that F = F\. We will show that F\ satisfies also the remaining 
conditions of Theorem 11.11 

By Corollary 4.8 , Fx = F* — G has the form (1.2) with completely mono- 
tone G\{x) = G{x), and G\ is the Laplace transform of a finite measure 7^ 
on (0, 00). By Lemma 4.9, Fx is the distributional eigenfunction of A(o^oo) - 

If ip is unbounded, then Fx = F vanishes continuously at (Lemma 4.11 ), 
and lima;_s.oo(-F\(a;) — sin(Ax -|- ^x)) = (because Gx = C'jx tends to 
at infinity). Furthermore, the subordinate Brownian motion Xt has 
the strong Feller property, and is a regular boundary point of (0, 00) 
(see Preliminaries) 



4.2 



HO,oo)-rA 



ator on L'^ 
cases Fx is an L°° eigenfunction of 
Note that (see [25 J ) 

ft 



applies: Fa G "^(^(0,00) 
is a b( 
A by Proposition 



and 



Hence, Lemma 

: — ^/'(A^)Fa. When ip is bounded, then ^(0,00) is a bounded oper- 
'((0,00)), and Ai 



V^(A2)F 
(0,00) • 



4.4 



In both 



■°°)f(x) - F{x) = [ '~U(o,oo)F(x)t/s = -H^') [ P^''°^^f{x)ds 
Jo Jo 

for X G (0, 00). For a fixed x, this integral equation is easily solved, and we 



obtain Pj^^'"^^ F{x) = e ^^^^^^F{x), that is. Fx is the eigenfunction of P^ 



(0,00) 
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B y Pr oposition 4.7, '&x = Argilj\(i\). Recall that tpxi^) = ^ O^J Proposi- 
tion k3^d)), and hence (by representation (2.23) for ■ip\) Iieip\{iX) > 1. It 



follows that < '(9a < 7r/2, as desired. It remains to prove formula (1.4). 

Suppose that (the holomorphic extension of) extends to a continuous 
function ip^iC,) (not to be confused with ip'^i^)) in the region Im^ > 0, and 
^+{-^) ^ V^(A2) for all ^ > 0. Recall that CGx = g (w ith g given by ( |47II| ), 
and Gx = C'jx- Hence, by the remark following (2.24), 7a = 7r~^mo, where 
rfio is the representation measure in (2.24) for the Stieltjes function g. To 
prove formula (1.4), it suffices to apply Proposition 2.18 , as we now describe. 

We have 



9{o = no-fio = no 



By Lemma 3^, when Im^ > and ^ ^ iX, 



9(0 = no- 



no - 



v^i(-0(A2+e) " A2V;i(-o(^(A2)-v^(-e)) 



Hence, g extends to a continuous function g~^ in the region Im^ > 0, and 
for ^ > 0, 



9^-0 = n-o- 



cx^iX') 



A2^1(0(^(A2)-^+(-a) 



Note that /*(— ^) and ip\{C,) are real for ^ > 0. By Proposition 
{—lmg~^{—C,))dC, is the measure rfiQ in the representation (2.24) o: 
Stieltjes function g. Since 7a = vr'^mo, Cx 
^(J^(X^) = 7/^(A2)/(AV(A^)), we obtain 



2.18 



the 

v/AV(A2)/^/'(A2) and 



lmg+i-0 



n 



1 VV''(A2)^(A2 



^ ^liO 



Im 



1 



Im 



^(A2)_^+(_^2) 

A^'(A2) 



d^ 



^(A2)-^+(-a 



dC 



An application of (4.9) completes the proof. 



□ 



Remark 4.12. Using (4.9), formulae of Theorem |1.1| can be written in a 
more concise form. We record them for easier reference: 



(4.13) 
(4.14) 

(4.15) 



CFxiO 



X 



X' 



A^-e 
,t 



log -7Y^ d( 
^a(A2) 



Arg^i(iA); 



,,,, 1 ^^^ , AV^XA^) 

^^^'^^^^^I(^'"^(A^)-^-(-e) 



d^. 



□ 
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Remark 4.13. The formula (1.4) for 7a remains true for general CBF ip, 
as long as it is understood in the distributional sense. More precisely, it can 
be proved that 7a is the limit of 



as e — > 0+. We omit the details. □ 

Remark 4.14. In order to extend Theorem |1.1| to more general symmetric 
Levy processes, one needs two ingredients: formula (1.1 ) must define a func- 
tion Fa, and the convolution F * *S)* i?_ needs to be associative (that 
is, some regularity of F and R± is needed). This seems to be problematic, 
even for general subordinate Brownian motions. □ 

4.4. Further properties of eigenfunctions. We begin with a technical 
result, useful for the computation and estimates of dx. 



Proposition 4.15. We have 

1 /"^ — logz J vr 

Proof. By series expansion and integration by parts. 



(4.16) - / -^^dz 



1 

2n 



/ ^^dz = y z'-{-\ogz)dz = y^— 
(2n + 1)2 ^ 1^1^' 1^ (2kY ^ 4^ P " 4 Y' 

n=0 ^ ^ k=l k=l ^ ' k=l 

as desired. □ 



The following formula for i^a is often more convenient than (1.3). 
Proposition 4.16. For A > 0, 

(4.17) ..^ir^iog*^.*.. 



Proof. By substituting = in (1.3), we obtain 

I r I , AV(A')(i-e') 

A substitution ^ = z in the integral over (0, 1) and ^ = I/2; in the integral 
over (1, cxd) yields 

^A = ^ ^ log dz 
^ n Jo 1-z^ ^ (^(AVz2) _ ^(A2))22 



Proposition 4.17. The function i^a is differentiahle in X > 0, and (4.17) 
can be differentiated in A under the integral sign. 
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Proof. Let h{s,z) = {\og{'ijj{s) — ip{sz'^)) — \og{^p{s/z'^) — ^{s)) — 2 log 2; for 
z e (0,1) and s > 0. Then h{X'^,z) = \og{ipx{X'^z^)/^x{Xyz^)). Therefore, 
it suffices to prove that {1 — z'^)~^\{d/ds)h{s, z)\ is integrable over rectangles 
z G (0, 1), s G (si, S2), where < Si < S2- 
We have 

dh{s,z) 1 siIj'{s) — sz^iIj^sz^) 1 sz~'^il!'{s/z'^) — sijj'{s) 
ds s ipi^s) — ip{sz^) s ip{s/z'^) — ipi^s) 

Using twice Cauchy's mean value theorem (for two functions til)\t) and 
ip{t)), we obtain that for some ti G [sz"^, s] and ^2 ^ [s, s/z"^], 

dhjs, z) ^ 1 ^'{t^) + t^r{t^) 1 ^'{t^) + t2r{t2) 

ds s ij'{ti) s ^'(ta) 

Since t^jj" (t) / ijj' (t) is bounded (by Proposition 2.21[ b)), the above expression 



is uniformly bounded in s G (si,S2) and z G (0,1). Furthermore, when 
s G (si, S2) and z G (1/2, 1), by the usual mean value theorem. 



dh{s^ z) 



ds 



< sup 

s te[sz2,s/z2] 



d fit) + tip" (t) 



dt 



< 



1 



sup 

iG[si/4,4s2] 



ijj'it) 

d ip'{t) + tip"{t) 



dt 



iP'it) 



< C{l~z' 



for some C depending on Si and S2- The proof is complete. 
Proposition 4.18. We have 



□ 



(4.18) 



< sup 



Proof. Note that (d/ds) log ipx^Xh) = \'^ij\{Xh)/iJx{\'^s). Hence, by ^J7\ 

"1 



It follows that 



1 



1 



ds dz. 



mo 

sup , /.X 

5eR V'a(0 



1-^2 



ds dz. 



Finally, by (4.16) 



1 



■'/^'l , ^ , 4 f' -logz ^ n 

ds I dz = — I ^ dz 

\ — z'^ 2 



TT 



□ 



Corollary 4.19. If = ^ for some A > 0, then ip is linear, Xf is a 
Brownian motion, and i)x = for all X > 0. 

Proof. Since ipx is increasing, for z G (0, 1) we have ipxi^X^ I z'^) > ipxi^X^z"^), 
and it is easily checked that equality holds if and only if ip is linear. Hence, 
by (4.17), dx > 0, and if ip is not linear, then dx> ^. □ 



Now we turn to the properties of Fx 
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(4.19) 



lim CFxiO 

f-s>0+ 



A > 0. 



Proof. The result follows from 4'\{0) = 1 (see Proposition 3.4 |^d)), CFx{^) 
CAV'UO/lA' + e'), and(|g. 

Lemma 4.21. For all X,x > 0, < Gx{x) < sin-i^^- Furthermore, 

PCX) 

(4.20) / Gxix)dx 



□ 



A 



COS-t?) 



'AV(A 



A > 0. 



Proof. Clearly, < Gx{x) < G'a(O), and by Lemma 4.11, we have < 
lima;_!.o+ Fx{x) = sin'dx — G'a(O). This proves the first statement. Further- 
more, CGx{0 = 9{i) = f*{0 ' fiO = - ^(0- By (pi), KM, 



(4.19) and Proposition 4.7 



POO 

/ Gxix)dx = lim CGxiO = /*(0) - lim CFxiO 



cxReijliiX) /V''(A2) cos^^a N^'{>^^) 



A2 

Proposition 4.22. We have 

|A + ei 



V'(A2) 



A 



^(A^ 



□ 



(4.21) 



A > 0, Re^ > 0. 



iA2+ei' 

Proof Since £Fa(0 = ca^1(0/(A^ + ^^), by (^ and (|48|) we have 

A 



l^i^A(OI < 



iA2+e 



1 



By Proposition |2^b), < 2^'(A2), and ( |42l] ) follows. □ 

Lemma 4.23. Suppose that there are G,a > such that > C^"^^ for 
^ large enough. Then (of. (4.21 )j 

(4.22) icFxiai < A /sup ':, .J.:::.: i + 



.r./.M / A2-2a^'(A2) A 



1-a 



whenever A > 0, and Re^ > 0. 

Proof. Since is concave and increasing, for A,.^ > we have 

^(A^) A^-e ^^(A^) / 1 1 
< — — — max 



MO 



A2 ^(A2)-V;(0- A 



Let Cx be the supremum in (4.22) (finite by the assumption), so that 
ij'{X'^)/f{0 < Cxi^/X'^y-'' for all ^ > A^. Then 



MO < 



AV(A2) 



1 + 
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By Proposition |3.5[ we have 
(4.23) 



l-a 



A > 0, Re^ > 0. 



The lemma follows from £Fa(0 = ca^1(0/(^^ + A^) and (4.8). 



□ 



Lemma 4.24. Suppose that there are C,a > such that ip'{C,) > C.^"^^ for 
sufficiently large ^. Then for each A > 0, the function Fx is Holder con- 
tinuous with any exponent less than a. More precisely, there is an absolute 
constant c such that 



(4.24) 



A2-2"^'(A2) 



a — e 



\Xi - X2| = 

for all £ G (0, a), A > and Xi, X2 G R 



sup 



OA C'-'-^'iC) 



Proof. We decompose CFx = CaV^^IO/IA^ + ^^) into the difference of /* and 
^ in a similar way to the decomposition CFx = /* — (? in the derivation 
of the formula for Fx, but this time /* and g have a better decay rate at 
infinity. When Re^ > 0, define 



cx 



A^l(zA) 



A^I(-^A) 



2zX \^(e-a)(e-2A + A) (e + a)(e + ^A + A) 



2iX\ - iX){^ - iX + X) 

(e + 2A)(e + a + A) 



Note that A(^ =F iX) ^ (^ =F iA + A) ^ is the Laplace transform of the function 
(1 - e-^^)e^*^^l(o,oo)(a;), and cx^i{±iX) = Xe^'^\ Therefore, /* is the 
Laplace transform of 

F*{x) = (1 - e-^") sin(Ax + ^9a)1(o,oo)(x). 

Note that F* is Lipschitz continuous on R. In fact, for x > we have 



\{F*y{x)\ = |Ae-^^sin(Ax + ^9A) + A(l -e-^^)cos(Ax + t9 
< v^A2e-2^^ + A2(l-e-^^)2 < A, 



so that the Lipschitz constant of F* is not greater than A. 

The function g extends continuously to zR. Denote the summands 
in the parentheses in the definition of g by ^+ and so that g{C,) = 
(ca/(2zA))(^+(0 -g-iO)- We first consider g+{i^) for ^ G [A/2,2A]. We 
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have 



< 



< 



{t^ - iX){i^ - iX + X) 

x\^i{^x) - ^im + \^ - xw^im 

ie-A|v/(e-A)2 + A2 

1 



A sup Mntc)\ + \i^{m\ 



y(e-A)2 + A2 

when A > ^, t hen we understand that [A,^] denotes the interval A]. 



By Proposition 2.21 c,d), |(^i)'(^C)l < 2^i(C)/C and |^1(2C)| < 2^I(C) for 



C G R; the latter inequality is used frequently in the remainder of the proof. 

6V^1(2A) 



Using also \/{^ — A)^ + A^ > A, we obtain 



\9+{i0\ < T A sup 



2^1(0 



A \ C6[A/2,2A] C 



A 



By (l423| , caV'1(2A) < S^^^Av^ < 3X^/C^, where Cx is the supremum 
in (itm Hence, 



A > 0, ^ G [A/2,2A]. 



When ^ G [-2 A, A/2], we simply have 1^ - A| > A/2, - A)2 + A^ > A, 
so that 



A|^l(zA)| 



ie-A|v/(e-A)2 + A2 |^-A| 



I 1^1(^01 ^4^1(A)+4^I(2A) 



A 



Using ( |423| ), we obtain that cx^liX) < 2i-"Av^ < 2Xy^. Together 
with c\iIj\{2X) < 3Xy/C\, this gives 

\cx~g4^0\ < 20v^, A > 0, e G [-2A,A/2]. 

Similar estimates hold for ca^_. It follows that 



|~^^^)| < c,i\~g+m + \9-m) < 20v^ 



2A 



A 



A > 0, ^ G [-2A,2A]. 



The estimate for ^ G R \ [-2A, 2A] is obtained using \g{iO\ < \CFx{iO\ + 
(here CFx{iC,) denotes the continuous boundary limit). By the in- 
equalities 1^ - A| > 1^1/2, - A)2 + A2 > 1^1/2, we have 



+ 



IV^1(-^A)| 



2 Vl^-A|\/(e-A)2 + A2 |e + A|v/(e + A)2 + A2 



^ 8cA^l(A) ^ 16AvC\ 
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Hence, using also (4.22) and the inequality — > (|^| + A)^/3, for 
^ G R \ [-2A, 2A] we obtain 



1^(^01 < 1^(^01 + ir(^oi< 

3V^ 



Av/gI(l + |e|/A) 



1+a 



This way we have estimated ^ ^ R- The function g is the 

Laplace transform of G, and G is the inverse Fourier transform of g{—iC,) 
e R). Fix e e (0, a). We have 



\G{xi) - G{x2)\ 



27r 

1 

< -TT sup 

2vr 



Since |e '^^^ — e < min(2, |^| |xi — X2I), the supremum is at most 

2^~'^\xi — X2f < 2|xi — X2f- Furthermore, 

+ 



20a/Ca 



1 + £ « — £ 



a — e 



a — e 



In particular, G is Holder continuous of order £, with Holder constant at 

□ 



most (200/71) A^a/C*a/(« — s)- The lemma is proved. 
Corollary 4.25. Suppose that for some e > 0, 



(4.25) limsup^^^ < 1 



hmsup <l-e. 



(see (l.ll)j. T/ien i/iere is C > Q (depending on e) such that 

\Fx{xi) - Fx{x2)\ < CX'\xi - X2|^ A > 0, xi,X2 G R, 

and 



\CFxm<C 



A 



A2+e 



l-e 



A > 0, Re^ > 0. 



Note that in general, < 2^/''(^), see Proposition |2.21[ b). 
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Proof. Choose a > e small enough, so that (4.25) holds with e replaced by 
a. Let h{C) = C^"°^'(C)- For ^ > A > we have 



exp 



exp 




dc 



X- HO 

\rio\ - 



dC 



The integrand on the right-hand side is positive for ( G (0, Ci) U (C2, 00) 
for some Ci,C2 > 0. Hence the above expression is bounded uniformly in 
A and ^, < A < ^. The result follows by Lemmas 4.23 and 4.24 □ 



Recall that a function / is said to be regularly varying of order a > at 00 
or at 0, if for all k > 0, lima;^oo f{kx)/f{x) = or lima.^o+ f {kx) / f {%) = 
k", respectively. 

Corollary 4.26. Suppose that ip is regularly varying of order Ooo > at 00 
and regularly varying of order ao > at 0. Then for all e < min(ao,aoo); 
the hypothesis of Corollary 4 -25] holds true. 

Proof. Let e < a < min(ao,aoo)- Since ip is a regularly varying Bernstein 
function, also tp' and ijj" are regularly varying at and 00 (see |6]). Hence, if 
we define h{() = C^~"'0'(C); then h is regularly varying of order ao — a > 
at zero, and regularly varying of order — a > at 00, and a similar 
statement with orders decreased by 1 is true for h'. It follows that the 
function (h'{()/h{() (continuous on (0, 00)) has positive limits at and 00, 
and we can now repeat the proof of Corollary |4.25[ □ 



When ip is regularly varying at infinity, then the behavior of Fx near 
for a fixed A > has a very simple description. We write f{x) ~ g{x) as 
X — )■ O"*" or X — > 00 when f{x)/g{x) converges to 1 as x — > 0^ or x — )■ 00, 
respectively. 

Lemma 4.27. Ifip is unbounded in (0, 00) and is regularly varying at 
infinity of order a G [0, 1], then 



(4.26) 



(4.27) 



Fx(x) 



'AV(A2 



r(l + a) 



'AV(A2) 

'ijj{x~'^) 



as ^ ^ 00, 



as X — >■ 0" 



Proof. Consider an auxiliary function: 

me) 



2/-2\ 



e,C>o. 



Since ip\ is a complete Bernstein function, both ipxlO and C,/'ip\{C,) are 
increasing (Proposition 2.17 a)). It follows that for ^ > 0, ( > 1, we have 
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Me)<Mee)<eMe)- Hence 



(4.28) 



<Me,C)< 



^(^2)^2-2. ^2 -■----^(p)^-2a ^2 

Note that ip{i'^)i>\{i'^) / i'^ has a positive hmit ip{X'^)/X'^ as ^ — oo. Hence, 
it is bounded above and below by positive constants when ^ G (l,oo). It 
follows that for some C > (depending on A), 



iogMe,C)l<C(i + |iogC|) 



e>i, C>i- 



For ( G (0,1], the inequalities in (4.28) are reversed, and we obtain the 
same bound for \ \ogh{C,,C)\- Since ip is regularly varying of order a and 



unbounded, 



lim h{C,, C) = lim -yr^ 



2a2 



lim 



xym' - 1 



for all C > 0. Hence, by dominated convergence. 



lim 



l0gV^A(e\'^ 



lim 



1 + e 

'log/^(^,C) 



TT 



1 + C 



2^°s A2^(e) 



Since the integral of (logC)/(l + over (0, oo) is zero, we have 



1 + C 



^(A2)C2-2a^2 

2^°^ AV(a 



1 ^ 

2 ^A2^(e) 



These two formulae yield that (see (3.1)) 



exp 



i + C^ 



A2^(e) 



as ^ — !■ oo. 



Since CFx = cxipliO / i^"^ + A^), this (and (|48|)) gives 



I jj'jX^) A' 



/^(A2)e 



A2^(e 



as ,^ — 7- oo. 



and (4.26) follows. Since Fa is increasing on the interval (0,7r/2 — -^x) 
(see (1.2)), formula (4.27) follows by Karamata's Tauberian theorem and 
monotone density theorem (see [6j). □ 

5. Spectral representation in half-line 

In this section we study the L'^{{0, oo)) properties of the operators pI^^'°°'' 
and prove Theorems 1.3, 1.7 and 1.8 Let us define the operator H* by the 
formula 



(5.1) 



fiX)Fx{x)dX, 



X > 0, 
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for / e Cc((0, oo)). Here Fx{x) = sin(Aa: + -dx) — G\{x) is given by Theo- 
rem [lT] At least formally, 11* is the adjoint of 11 defined in ( 1.5[ ). 

In this section, we write ||/||2 for ||/||l2((o,oo)), and {f,g) for (/, 5')l2((o,oo))- 

Lemma 5.1. The operator U* extends to a bounded operator on L'^{{0, oo)). 

Proof. We follow the proof of Theorem 3 in [49j. Let / G Cc((0, oo)). Define 

poo _ 

n*/(x) = / f{\)sm{Xx + ^x)dX = lmJ^f{x), x > 0, 

Jo 

where /(A) = e^'^^fiX). Clearly, ||n*/||2 < WJ^/WlHr) = V^ll/lh- Let 

POO 

n;/(a;)= / fiX)Gxix)dX, x > 0, 

Jo 

so that n* = - n^. We have 

) POO POO POO 

iU*J{x))^dx< / / / \fiX,)\\fiX2)\GxAx)Gx,ix)dxdX^dX2. 
Jo Jo Jo 



By Lemma [4.21 , when < Ai < A2, 

Gxi{x)Gx2ix)dx < sini^Ai / Gx2ix)dx 



By symmetry, for all Ai, A2 > 0, 

Gxi{x)Gx2{x)dx < 



max(Ai, A2) 

Hence, using Hardy-Hilbert's inequality (see e.g. [33j) in the last step. 



1=1 m{x)fdx<l I ""-^^^dX^dX^KAWnl 



The lemma follows. □ 

The extension of H* to L^((0,oo)) is denoted by the same symbol. By 
Lemma 5.1, we see that H* is indeed the adjoint of H defined in ( |1.5 ). 



Lemma 5.2. We have {Wf.Wg) = (7r/2)(/,c/) for f,g e L^{{0,oo)). Fur- 
thermore, for f E L^((0,oo)) such that e~^^^'^ ^■^ integrable in X > 0, 
we have 

POO 

(5.2) '°"^n7(x) = / e-*^(^'V(A)FA(a;)rfA, t,x>0. 

Jo 

Proof. Again the argument follows the proof of Theorem 3 in [49j. Let 
pf''°°\x,dy) be the kernel of pj:^'°°^ (i.e. the P^; distribution of the killed 
process Xj[°'°°''). If / G Cc((0, 00)), then pf''°°\x, dy) f {X) Fx{y) is integrable 
my,X > 0. Theorem 1.1 and Fubini yield ( |5.2 ) in this case. The general 
case / G -^'^((O, 00)) follows by approximation. 
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Let f,g e C,((0,oo)), k e Z, and define MX) = e-^*^(^')/(A), 
g^(X) = e-'=*^'(^')^(A). From 1^ it follows that '°°^n*/fc = Tl*fk+i 
and pf '°"^n*(7fc = n* gk+i- The operators p^(°'°°) are self-adjoint, so that 

n*f_,,Pi''^^U*g) = {Il*f^,,U*g,). 



By induction, we have {Il*f,Il*g) = {U* f_k,Il*gk) for k > 0. In par- 
ticular, if supp/ C (0, Ao) and suppt? C (Ao,cxo), then {Il*f,Il*g) = 
l^^-ktiij{\) j_^^^kti){\) The right-hand side tends to zero as /c — > oo, 
so that n*/ and TVg are orthogonal in i^^((0, oo)). By approximation, this 
holds true for any f,g & L^((0,oo)) such that /(A) = for A > Aq and 
^(A) = for A < Aq. 

Define m{E) = llnn^^lH for Borel E C (0,cx)). Clearly, < m{E) < 
C\\1e\\1 = C\E\. If Pi C (0, Ao) and P2 C (Aq, 00), then 

m(PiUP2) = \\U*lEAl + WlE,\\l + 2{U*lE„Il*lE,)=m{Ei)+m{E2). 

Finally, suppose that P = IJ^i-^w where Pi C P2 C ... and |P| < 00. 
Since 1e„ converges to 1e in L^((0,cxd)), also n*!^;^ converges to n*!^; 
in L^((0, 00)), and so m{E) = lim„_>oo "^(-^'n)- It follows that m is an 
absolutely continuous measure on (0, 00), and by approximation, 



/■oo 

{n*f,U*g)= / fi\)gi\)mid\) 
Jo 



for any f,g & L^((0, 00)). The lemma will be proved if we show that 
m(P) = (vr/2)|P|. 

Fix Ao > and define fs = (1/v^)1[Ao,Ao+5] ioi 6 > 0. Clearly, Wfsh = 1, 
and m([Ao, Ao + S])/6 = ||n*/5||2. Furthermore, VSU*fs = gi + gi — 93, 
where 

pXq+S 

gi{x) = / sin(Aa; + 9xo)dX, 

pXq+S 

g2{x) = / (sin(Ax + 9x) - sin(Ax + 9xo))dX, 

gz{x) = / Gxix)dX. 
We study the behavior of gi, g2 and g^ as S ^ 0^. 



By Lemma 4.21, we have Gx{x) < 1 and Gx{x)dx < 1/A. Hence, by 
Fubini, 



^73||2 = / ( / Gx{x)dX ] dx< I 6[ I Gx{x)dX ] dx < 



Xo+S \ 2 /-oo / rXo+S \ r2 



Ao' 



'Ao / Jo \JXo 

SO that ||5'3||2/v^ tends to zero as 5 — )■ 0^. 

The function (72 is the imaginary part of the Fourier transfor m of the func- 
tion /(A) = (exp(i7?A) - exp(i^9Ao))l[Ao,Ao+<5]- By Proposition 



4.17 



we have 
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[''^A— ''^Aol ^ Ao| for some C > (depending on Aq). Hence |/(A)| < C5, 
and so ||/||2 < C^^/s. it follows that Hc^alh/v^ < v^H/lb/v^ < v^C(5 
also tends to zero as 5 — )■ 0"*". 

Finally, for the function gi, by Fubini, 



igiix)ye-'^dx 

Ao+<5 Ao+5 oo 



Ao Ao 
Ao+<5 Ao+(5 oo 



J 2 sin(Aix + 6'ao) sin(A2a; + 9\^)e '^^dxdXidX^ 



Ao Ao 
Ao+<5 Ao+(5 



(cos((Ai — X2)x) — cos((Ai + X2)x + 26xo))e ^'^dxdXidX2 
£cos(2^?Ao)-(Ai + A2)sin(279Ao^ 



Ao Ao 



+ (Ai - A2)^ 



+ 



+ (Ai + A2) 



dXidX2. 



By taking the limit e — )■ O"*" and using dominated convergence for the second 
term of the integrand, we obtain that 



I l|2 



{gi{x)) dx 



Ao 



7r 



sin(2^9Ao) log 



Hence, using < log(l + s) < s (s > 0), 



lim 

5^0+ 



I l|2 



TT 

2 



lim 

6 2 5-^0+ 26 
We conclude that for all Aq > 0, 

m([Ao, Ao + 5]) 



sin(2t9 



Aoy 



Ao+<5 



Aq 



log 1 + 



Aq + ^ + A; 
An + A2 



dXo 



An + As 



ciAo. 



TT 

2 



lim 



6 



lim 

5^-0+ 



191 + 92 -93h 



IT 

2 ' 



that is, m{dX) = {n/2)dX. 



□ 



5.2 



2/7rH* is an 



Proof of Theorem | when H zs injective. By Lemma 
L^((0, 00)) isometry, and by the extra assumption, its adjoint A/2/7r H is 
injective. Hence, both operators are unitary on L^((0,oo)). The theorem 
follows now directly from Lemma 5.2 and (H*)^^ = (2/7r)n*. □ 



In [56j, Theorem 1.3 is proved in full generality. More precisely, it is 
shown that H is always injective. Here we are satisfied with the following 



observation: When Xt is a symmetric a-stable process, then, by Example 6.1 



below, Fx{x) = -Fx (A), and therefore H* = H. Since a self-adjoint isometry 
is necessarily unitary, we see (by Lemma 5.2 ) that H is injective. Hence, we 
have proved the following result. 

Corollary 5.3. When a G (0, 2] and Xt is the symm etric a-stable process 
(that is, = ^"^'^), then the hypothesis of Theorem 1.3 holds true. □ 
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In the following proofs of Theorems L7 and L8, we use the full statement 
of Theorem 1.3 In other words, we assume that 11 is injective (which is 
proved to be true in general in [56j). 



Proof of Theorem \l .'T\ Let / G Cc((0, oo)). By Theorem 1.3 
2n*(e-*'^(A^)n/(A))(x' ^ 



p. 



(0,oo) 



TT 



e-'^^'^')Fx{x)irf{X)dX 



for all t,x > 0. The function e ^'^'''^^^ F\{x)F\{y) f {y) is jointly integrable in 
I/, A > 0. By Fubini, 



,(0,oo) 



oo /2 
TT 



Since / G Cc((0, oo)) is arbitrary, the kernel of pj^^''^'^ has the form (1.6). □ 



Proof of Theorem 1.8 under additional assumption (1.11). By monotone 
convergence. 



P.(r(o.oo) >t) = P^'^^h 



(0,oo) 



X) 



lim P, 



(0,oo) 



eAx] 



where e^{x) 



e '^^l(r)^\{x). We have 



L(0,oo)l 

ne,(A) 



e-'''Fx{x)dx = CFx{e). 



By Theorem 1.3 
(5.3) 1 



,(0,oo) 



ee{x) 



/■oo 

/ e-'^^'"^F,ix)CFxie)dX. 
Jo 



2 

71 



If we can use dominated convergence, then, by Proposition 4.20 

2 



Pa-(7-(0,oo) > t) 



-m^')F^{x) 



limCFxie) 1 d\ 



rfA, 



as desired. Hence, it remains to show that the integrand in ( 5.3^ is domi- 
nated by an integrable function when £ — )■ 0+. 

By the assumption (1.11), th ere a re C,a > such that ipiC,) > C^" for 
^ > 1 (see the proof of Corollary 4.25). Hence, e~^'^^^ -* is integrable in A > 0. 
By Proposition 4.22, we have CFx{e) < 2/A, and by Corollary 4.25 (for 
Xi = and X2 = x), \Fx{x)\ < C min(l, A"/^) for some C > (depending on 
x). Hence, the integrand in (5.3) is bounded by 2Ce~^^^^ min(l, A°/^)/A, 
which is integrable in A > 0. Formula (1.9) is proved. 

By the same argument, the integrand in (1.10) is integrable jointly in 
A > and t > t^. Formula (1.10) follows by integrating both sides of it in t 
over (to, C)o), using Fubini and comparing the result with (1.9). □ 
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Example 6.1. Let tp{^) = ^"/^ where a G (0,2). Then Zt is the (a/2)- 
stable subordinator, and Xf = Bz^ is the symmetric a-stable Levy process. 
By ( |4l7| and ( |4l6| ), 



log 



dz 



2 — a /"^— log2 , (2 — a)7r 
dz = 



(Note that the phase-shift of (2 — a)7r/8 was conjectured, in somewhat 
different setting, in |73j; see |54j for further discussion.) By Theorem |L1| 
the eigenfunctions of pj^^'^^"^ and ^(o.oo) are given by the formula 



POO 

Fxix) = sin(Ax + (2 - a)n/8) - / ^xiOe'^'^d^, 

Jo 



X > 0, 



where (see (4.15) and (3.1)) 



7a(0 



/2aA 



a-l 



sin(a7r/2) 



2tt A2" + _ 2A"^° cos(a7r /2) 



X 



X exp 



log 



i + C i-ec°A" 

Clearly 7(5) = A7a(As) does not depend on A, and finally. 



dc 



where 



Fx{x) = F{Xx) = sin(Ax + (2 - a)n/8) 



7(s) 



7(sje 



-\sx 



ds. 



2a sin(a7r/2) 



27T 



X exp 



1 + - 2s° cos(a7r/2) 

— ■ log :; — dC 



1 + C 



By Lemma 



4.27 



we have Fx{x) ~ (v/^r(a/2))-i(Ax)"/2 as x ^ 0+. 
By Theorem 1.3 (or Corollary 5.3), the functions Fa yield a generalizec 

and 



eigenfunction expansion of ^(0.00) and 

2 

IT 



{0,00) 



and by Theorems 



1.7 



(0,00) . 



Px(T(0,oo) > t) 



-tX°' 



F{Xx)F{\y)d\, 





/2a 



TT 



00 g-tA° 



F{Xx) 



A 



dX, 



For a 



1, these results were obtained in 



t,x,y> 0, 



t, a; > 0. 



and the formulae for 



(0,00), 



x,y) and Px(t"(o,oo) > t) can be significantly simplified in this case. 



Pt 



Example 6.2. If ip{^) = \/m^ + ^ — m for some m > 0, then is the 
relativistic subordinator with mass m, and the corresponding Levy process 
is the relativistic 1-stable process with mass m. Its infinitesimal operator A 
is the (quasi-)relativistic Hamiltonian of a free particle with mass m, and 
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^(0 oo) is the Hamiltonian corresponding to the semi-infinite potential well. 
By 'dm, 



TC 



log 



(^y/m^ + Xyz^ - Vm2 + A2 j z^ 



dz. 



which increases with A from as A — ?■ 0^ to tt/S as A — ?■ oo. The eigenfunc- 

sin(Ax + {}x) — Gx{x), where G\ is 



tions of P4 



(0,00) 



have the form F- 



the Laplace transform of an explicit nonnegative function. By Lemma [4.27 
TT as X — )■ 0"''. See I 

example. 



for a detailed discussion of this 



Example 6.3. Let ipi^ = + ^ (o,2), /3 > 0. Then Zt is the 

(a/2)-stable subordinator with drift, and the corresponding Levy process 
Xt is a mixture of the Brownian motion and the symmetric a-stable Levy 
When a = 1, then, by (|4.17D, 



process. 



1 



TC 



TC 



, (A + /3A2) - (AC + PX^e) 

log . R^2^ . ....... dC 



1 



log 



, (A + /3A2))C2 
l + f3X + f3X( 



l-C "C + /3A + /3AC 



dC 



This decreases with A from vr/S as A — )■ 0^ to as A — )■ 00, and can be 
written explicitly in terms of the dilogarithm function Li2. For general a, the 
expression for 'dx is more complicate d, and it can be proved that -t^x decreases 



from (2 — a)Tc/8 to 0. By Lemma 
X 0+. 



4.27 



Fx{x) ~ + (a/2)A"-2 Xx as 



Example 6.4. Let ip{C,) = log(l + ^), so that Zt is the gamma subordinator. 



The subordinate process Xt is called variance gamma process. By (4.17) 



^9) 



TC 



1 



log 



log(l + A2) - log(l + X^z^) 

(log(l + A2/^2) _ log(l + A2))z2 



dz . 



It can be pr oved t hat 'dx increases with A from as A — ?■ 0+ to 7r/4 as A — ?■ 00. 

4.27| the eigenfunctions Fa satisfy Fx{x) ~ Cx/ a/| loga;| as 
with Cx = X/ \J2 + 2A2. In particular. Fx is not Holder continuous. 



By Lemma 
X 

Note that the assumptions of Theorem L7 are satisfied only for t > 1/2, but 



Theorem 1.8 applies for all t > 0. On the other hand, the more restrictive 



condition (1.11) is not met, so this case is not covered by the proof of 



Theorem 1.8 given above. The general argument of [56] needs to be used. 



Example 6.5. Let '?/'(0 = log(l + log(l + ^)). It can be verified that in this 
case dx is greater than 7r/4 for some A, e.g. "d^ ^ 0.2877r. This proves that it 
is not true in general that ^x < 7i"/4, even if ip is unbounde d. Fu rthermore, 
g-#(A2) integrable in A > for any t > 0, so Theorem 



1.7 



is of no use 



here. Theorem 1.8 applies for t > 1/2, but the additional assumption (1.8) 
is not satisfied. 
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Example 6.6. Let ipiO = ^/(^ + 0- this case the subordinator Zt is 
the compound Poisson process with exponential jump distribution, and the 
jumps of Xt have Laplace distribution with density (l/2)e~'^'. Note that 
l/('?/'(A^) — ipiO) is holomorphic in G C \ {A}, and hence 7a vanishes. It 
follows that Fx = sin(Ax + 'i?A)l(o.oo)(2;)- By ( 1.3[ ) and a contour integration, 



1 f°° A 1 + 

= — ^ log rfz = arctan A; 

27r A^ - 1 + A^ 

we omit the details. Hence, F\{x) = sin(Aa; + arctan A)l(o,oo)(a;) satisfies 
AFx{x) = —ip{X'^)Fx{x) for all x G (0, oo). This can be alternatively proved 
by a direct calculation. 

It can be easily proved that 7a vanishes if and only if ipiO — (i-^- -^t 
is a Brownian motion) or = Ci^/{^ + C2) for some Ci,C2 > (which 
corresponds to the process Xt studied in this example, up to time and space 
scaling). 

Example 6.7. In general, 7a may have non-zero singular part. For example, 
consider = 5x / {x + 1) + x / {x + 5) and A = 1. Then = 8/3 = 

and hence 7a has an atom at 2 for A = 1. 

7. Application to systems of PDEs 



As it was the case in j49|, where the Cauchy process was studied, the 
eigenfunctions Fa are the boundary values of solutions of a certain system of 
PDEs, a spectral problem with spectral parameter in the boundary. In [49j, 
the formula for Fa (for the case of the Cauchy process) was found by solving 
the corresponding system of PDEs. Here the argument goes in the opposite 
direction, and the result for PDEs is an application of the formula for Fa. 

We need some background on Krein correspondence. For a complete 
Bernstein function ip, there is a unique string m (a non-decreasing, right- 
continuous unbounded function from [0, 00) to [0,oo]) which is the Krein 
representation of ip. For simplicity, we only consider the case when m(0) = 
0, m is continuously different iable on (0, 00), and m'{y) > for all y > 0. 
Denote a{y) = l/m'{y). Then, for each A > there exists a (unique) 
nonnegative, non-increasing continuous function gx{y) = g{X,y) such that 
a{y)g"{y) = Xgx{y) for y > 0, and gx{0) = 1. Furthermore, 

^(A) = -g'xiO), A > 0. 

For the details, we refer to or Chapter 14 in |66) . and the references 
therein. The Krein's correspondence between m and ip is a rather compli- 
cated concept, and there are only few examples of explicit pairs, see [66]. 
Let / G 5 and define u{x,y) using the Fourier transform in x, 

FMU) = 9ie,y)J'f{0, e e R, y > 0. 

It is easy to check that u is the solution of the problem 

+ a{y) ^ ] u{x,y) = xeR,y>0; 



dx"^ dy 



M(a;,0) = /(x), a; G R. 
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Furthermore, 

\dy ) y^o+ y 

so that {d/dy)u{x, 0) = Af{x), where A is the operator with Fourier symbol 
Theorems imphes the following result. 

Theorem 7.1. Let a{y) be a positive continuous function on (0,oo) such 
that 1/a is integrable over (0,1), but not over (l,oo). For A > 0, let 
gx{y) = g{^,y) be the nonnegative, non-increasing continuous solution of 
0'i.y)9\i.y) = >^9\iy) (y > O), satisfying gxiO) = 1. Define ^p{X) = g'^{0), and 
let Fx be the function given in Theorem Finally, let 

(7.1) ux{x,y) = J^-\g{e,y)J'Fx{0), xeR,y>0. 

Then U\ is continuous and bounded in x & H, y > 0, twice differentiable in 
a; G R, y > 0, and it satisfies 

+ a{y) 7— ma(x, y) = x G R, y > 0; 



(7.2) MA(a;,0) = 0, x<0 

-^UA(a;,0) = -^(A2)ua(x,0), a; > 
dy 



The Fourier transform and inverse Fourier transform in (7.1) needs to be 
understood in the distributional sense. 

Sketch of the proof. Using the properties of g{X, y), it is easy to check that 
Ma is a weak solution of (7.2). Since gi^"^ ,y)J^F\{^) (for a fixed ?/ > 0) is a 



sum of an -^^^(R) function and Cy/ (^ — A) — Cy/ (^ + A) (in a similar way as 
in Section |4]), u{x,y) is a bounded function. Continuity of ux and its first 
and second partial derivatives in x G R, y > follows by elliptic regularity 
theorem. Continuity of Ux on the boundary is proved using the properties 
of g{X,y) and Fx{x). □ 

We provide two examples. If a{y) = l/{l + 2ay) for a constant a > 0, then 
we find that gx{y) = (l+2a?/)^(^^+^-")/(2a) andV'(A) = y/X+~a?-a. Hence 
ip is the characteristic exponent of the relativistic subordinator, described 



in Example 6.2 For a related example, see Section 2.7 in 

Consider now a G (0, 2) and a{y) = a^c^?/^"^/", where Cq, = 2~"r(l — 
a/2)/r(l + a/2). We have gx{y) = C^{c^^yfl^K^/2{{c^V\yf'''), 
where Ka/2 is the modified Bessel function of the second kind and Ca = 
Q;2~"/2/r(l + a/2). Moreover, il){\) = A"/^, so that V is the characteris- 



tic exponent of the (Q;/2)-stable subordinator, considered in Example 6.1 
When a = 1, a{y) = 1. 



The problem (7.2) with a{y) = 1 was studied in |49] to find the eigen- 



functions of p^^"'^ Earlier, a similar relation for more general domains 
(also in higher dimensions) was applied e.g. in |2l [3l [10], [70] , and for general 
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symmetric stable processes e.g. in [T9| l20l |2T| l59] . Related problems ap- 
pear frequently in hydrodynamics (the sloshing problem), see the references 
in i49j. 

The operator (P /dx^ + a{y)(P / dy^ is the generator of a two-dimensional 
diffusion on R x [0, oo), with reflecting boundary. By repeating the 
construction given in Example 3.1 in J42j, we can identify the subordinate 
Brownian motion Xt studied in the present article with the trace left on the 
horizontal axis by Yt. We remark that the connection between subordinate 
Brownian motion in R*^ and traces of diffusions in R'^^^ was applied e.g. 
in [T2| [13] to find formulae for the distribution of some Levy processes 
stopped at the time of first exit from the half-line or the interval. A related 
problem for the trace of a two-dimensional jump-type stable process was 
studied in [36] . 
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